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b n parameter in recurrence relation (see eq. (15)) 
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GAUSS QUADRATURE RULES INVOLVING SOME NONCLASSICAL 
WEIGHT FUNCTIONS 

Paul F. Byrd and David C. Galant 
Ames Research Center 


SUMMARY 


Gauss quadrature rules are given for numerically evaluating integrals on 
the interval [0, 1] for some weight functions of the form w(t) = tY (1 - t a )£, 
where a is a positive integer > 2 and 3 , y > -1. The nodes, weight coeffi- 
cients, and the parameters in the three-term recurrence relation for the 
sequence of orthogonal polynomials associated with w(t) are tabulated to 25 
significant figures for N-point rules with N = 2, 4, 6, 8, 12, 16, and 24. 
Numerical tables to 5S are also given for a coefficient in the expression of 
the error term. 


INTRODUCTION 


A wide variety of problems requiring numerical approximations to 
specific definite integrals arise very often in practice, usually as an 
intermediate step in a more extensive computation. There are many methods 
for approximating such integrals, but among all quadrature formulas of the 
form 


I 


J g (t)dt = J 


N 

w(t)f(t)dt - ^ W jN f(t jN ) 
j = i 


( 1 ) 


only the weighted gaussian types have all their parameters chosen to maximize 
performance. For these formulas, the nodes t-jN and the coefficients Wjn 
are uniquely determined so that all approximations to the given integral are 
exact for arbitrary polynomials f of the highest possible degree, 2N - 1 . 

If the function f is other than a polynomial of degree <2N - 1, the inte- 
gration formulas (1) have an error term given (e.g., see ref. 1) by 


where 


en 


k N 
(2N) ! 



a < x < b 



<f>N being the associated Nth-degree orthogonal polynomial. 


( 2 ) 

( 3 ) 


with 



The literature concerning weighted gaussian formulas is very large, both 
classical and recent (refs. 2, 3, and 4). Widespread use of high-speed com- 
puters, with emphasis on economizing the number of operations required in com- 
putations, has resulted in considerable renewed interest in the power of the 
gaussian technique; gauss quadrature rules are being tabulated and used with 
increasing frequency. Very accurate rules, for example, were recently tabu- 
lated in references 1 and 5 for cases associated with the classical orthogonal 
polynomials of Legendre, Laguerre, Chebyshev, Jacobi, and Hermite. Recent 
numerical tables to many significant figures have also been given for some 
nonclassical cases such as a weight function on [o, l] having logarithmic sin- 
gularities (ref. 6) and the weight function e" x2 on the interval [0, »] 

(ref. 7). 

Weighted gaussian integration formulas are excellent for large classes of 
functions arising in practice. They are particularly useful when a complete 
integrand g(t) is factorable into two functions f(t) and w(t), where f is 
well approximated by a polynomial, with the nonnegative weight function w 
having a singularity or a singular derivative at one or both ends of the 
interval. Since neither the error term (2) nor the sum in the approximations 
(1) directly involves the weight function w, once the quadrature rules are 
constructed and programmed on a computer, they can be repeatedly applied to 
integrals of the form (1) with different functions f(t) and the same w(t). 
There are several other advantages of gauss-type formulas: (a) since all the 

coefficients Wjjsj turn out to be positive, the rules have good roundoff prop- 
erties; (b) the sequence of formulas (1) for N = 1, 2, . . . will always 
converge to the value of the Riemann integral I whenever this integral 
exists; (c) the rules have more favorable error terms and generally achieve 
higher accuracy with the use of fewer nodes than Newton-Cotes methods. 

This report presents tabulated gaussian quadrature formulas for 

N 

I(a, g, y; f) = J w (a, g, y; t)f(t)dt - ^ ' W-j N f (t-j N ) (4) 

0 j = l 

where the weight function w is of the form 

w(a, g, y; t) = w(t) = tY(l - t a ) e (5) 


with a as a positive integer ^2; 3, y > -1; and w or its derivative 
singular at one or both ends of the interval. For the special case a = 1, 
the approximate relation (4) becomes the classical Gauss-Jacobi formula 

f tY(l - t)f*f(t)dt * (6) 

0 


2 


II 


Krylov, Lugin, and Janovic (ref. 8) recently gave these quadrature rules to 
eight significant figures for 8, y = -0.9 (0.1) 3.0, 8 < y and N = 1(1)8. 

We shall consider various cases where the parameters (a, 8, y) are 


a=2,3,4,6,8 ' 

8 = -3/4, -2/3, -1/2, 1/2 > 

y = -1/2, 0, 1/2 , 


( 7 ) 


One of these, for example, is the integral 

1(3, -1/2, 0; f) = f - f dt (8) 

J A - t3 

0 

which, when f(t) is simply a rational function, may be reduced by a substi- 
tution involving Jacobian elliptic functions to a combination of three stan- 
dard tabulated elliptic integrals (ref. 9). A similar statement applies to 
integrals I (a, 3, y; f) with regard to the other values a, 3, y given in 
equations (7) . 

If f is not a rational function (e.g., e~ ct ) , none of the integrals 
I (a, 3, y) are expressible in finite terms of known tabulated functions for 
the chosen parameters (ct, 3, y) , and the only recourse is numerical approxima- 
tion via some quadrature formula. Even when such integrals can be written as 
combinations of standard elliptic integrals etc., it is probably more effi- 
cient to use quadrature rules having the highest degree of algebraic precision 
to obtain answers. Integrals with our weight functions w(t) occur suffi- 
ciently frequently in practice to warrant the construction of such rules. 

Numerical tables for N-point gaussian rules {Wjn, tj^} with 
N = 2, 4, 6, 8, 12, 16, and 24 are given to 25 significant figures for 12 


special cases 

of I (a, 6, y; f) 

These 

are 

for 





f(3, -1/2, 0) 

(3, 1/2, 

0) 

(4, 

-1/2, 0) 

(4, 1/2, 0) ' 


(a, 8, y) = " 

(6, -1/2, 0) 

(6, 1/2, 

0) 

(8, 

-1/2, 0) 

(8, 1/2, 0) 

- (9) 


* (3, -1/2, -1/2) 

(3, 1/2, 

1/2) 

(2, 

-3/4, 0) 

(2, -2/3, 0) J 



The coefficients b n , g n in the three-term recursion relation (15) for the 
sequence of monic orthogonal polynomials corresponding to w(t) are also 
tabulated to 25S for N = 1(1)26. Since the error coefficients k^ in the 
remainder term (2) may be useful when a bound on the 2Nth derivative can be 
estimated, we give tabulations of them to 5S. 


3 


Little difficulty is entailed in the actual construction of gauss-type 
quadrature formulas relative to classical weight functions where the three- 
term recurrence relation for the associated orthogonal polynomials are known 
explicitly. This is not generally true for nonclassical w(t) such as ours 
where the sequence of related unknown orthogonal polynomials must be 
generated using the moments of w. 


GENERAL DISCUSSION OF CONSTRUCTION 


For a given nonnegative weight function w(t) on [0, 1], the problem of 
constructing the N-point gaussian quadrature 

N 

G N (f) = T: W iN f (t iN ) = J w(t)fCt) - E N (f) , (N = 1 , 2, . . .) 

j=l 0 (10) 

consists in determining the nodes tj^ and the weight coefficients Wjj^ (2N 
parameters in all) so that the error EN(f) = 0 whenever f is a polynomial 
of degree n < 2N - 1. There are essentially only two approaches to the 
determination of these unknown parameters, and they shall be discussed 
briefly. 

Before proceeding, however, we wish to emphasize in the strongest 
possible terms the numerical difficulties inherent in the construction of the 
quadrature rule (10) from the moments 


M 


k = 



t k w(t)dt , 


k = 0 (1) 2N - 1 


(ID 


Gautschi (ref. 4) has shown that an approximate upper bound to the relative, 
asymptotic condition number of this problem is 


„ (17 + 12/2) N 

“ 

64N 2 


( 12 ) 


This means that the computation is extremely ill-conditioned numerically, and 
that the number of significant figures lost in determining the rule (10) using 
the moments (11) is approximately 


Sn = logio C N = 1.53N - 2 logi o (8N) 


(13) 
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Thus, if one wants to have the nodes and weights of Gn accurate to k 
significant figures, then he must work with more than 


k + 1.53N - 2 log! o (8N) 


significant figures. For example, to obtain a 10-point quadrature rule with 
nodes and weights accurate to 10 significant digits (i.e., N = k = 10), one 
must work with more than 22 significant digits at a certain stage in the 
calculation of the nodes and weights. 


Algebraic Approach 


If the nodes tjN and the weight coefficients Wjn in equation (10) are 
chosen so that the formula is exact for all polynomials f(t) of degree 
n < 2N - 1, then it is sufficient to choose the formula to be exact for the 
monomials t^, k = 0(1) 2N - 1, so that one has the system 


N 

y; WjN tj N = M k , k = 0 (1) 2N - 1 (14) 

j=l 

of 2N nonlinear algebraic equations in 2N unknowns to solve, where 
are given by equation (11) . The drawback to this formulation of determining 
the quadrature rule from the moments is that any attempt at a direct solution 
requires the highly precise arithmetic indicated by equation (13) at every 
step of the calculations. For this reason, the algebraic approach is to be 
rejected except for very small N. 


Approach Via Orthogonal Polynomials 

Orthogonal polynomials play an important role in gaussian quadratures . 

In fact, it can easily be shown that if the nodal points t j jq in formula (10) 
are selected to be the roots of the Nth-degree polynomial cj)^(t) of the 
sequence of orthogonal polynomials (4> n (t)} associated with the given weight 
function w(t), then the formula is exact for any polynomial f(t) of degree 
< 2N - 1; the converse of this is also true (e.g., see ref. 10). 

The orthogonal system of polynomials associated with w(t) in (a, b) is 
unique and any three consecutive polynomials satisfy the recurrence relation 


4> n (t) = (t - bj^.^t) - g n V 2 (t) , n = 1 (1)N (15) 

with cj>_i(t) = 0, cj> 0 (t) = 1 , and 


5 



bn 



tw(t)<()^_ i (t)dt 
w (t)^_i (t)dt 


g 


n+i 



l 

w(t)cj> n 2 (t)dt 
w(t)*J_ i (t)dt 






J 


( 16 ) 


The usual approach in the construction of gaussian rules is based on the use 
of relation (15) and the moments to generate the sequence of polynomials 
{ (t ) } . After these are calculated, one solves, in some way, for the roots 
tjN of <j>jj(t), and then finds the weight coefficients WjN by other means. 
However, any procedure that relies on the moments is essentially equivalent to 
solving the basic algebraic system (14) and is highly unstable numerically. 

The advantage of the approach via orthogonal polynomials, as we shall 
see, is that the ill-conditioned nature of the problem is encountered in the 
determination of the three-term recurrence relation from the moments, but at 
no point in the determination of the nodes and weights of the quadrature rule. 
Since the calculation of b n , g n , n = 1(1)N, needs less than 4N 2 arithmetic 
operations, the amount of multiply precise arithmetic is small and manageable, 
thus making this approach very tractable in comparison to the algebraic 
approach . 


NUMERICAL METHODS 


The numerical construction of our gaussian rules proceeds in three steps: 

(1) Determine, to sufficient accuracy, the moments Mfc for k = 0(1)2N - 1; 

(2) using the moments, calculate the parameters bj^, g^, k = 1(1)N, in the 
three-term recurrence relation (15) for the first N orthogonal polynomials 
associated with the weight function w(t) on (0, 1); (3) determine the quad- 
rature nodes t-jN and the weight coefficients WjN for j = 1(1)N. The 
three steps will be discussed in that order. 


Determination of the Moments 
The moments for the weight function of interest are 


Mn (ct , 3, y) = M n = J t n w (t) dt = J t 


n+ ^(i - tV dt 


(17) 


H3 '0 

for n = 0(1)2N - 1. If we let z = t a , equation (17) becomes 
6 


Tlirs !" ■ ' “T"ll I '1 ■ 


■■ ■■■■ifMi ■■ ■■■■ i i -rr*n~r- i m hh ■■■ 



(a > 0; £5 , Y > -1) 


M n 



z (n+y+l)/(a-l) ^ _ z -j(l + 3)-i 


dz 


or 


^ -5- B 


n + y + 1 


1 + 6 


r(n + y + l/a)r(l + 6) 
aT[l + 6 + (n + y + 1/ct)] 


( 18 ) 


where B and r are the beta and gamma functions, respectively. Using the 
well-known recurrence relation for the gamma function, we have the recursion 
relation 


Mk = 


1 + [a (1 + 


M v 

K-a 

6)/(k 


Y - a + 1)] 


k = a (1) 2N - 1 


(19) 


so that all the moments can be generated from Mg , , . . . , ^a-1 


Calculation of Parameters in the Three-Term Recurrence Relation 


There are several methods of calculating the three-term recurrence 
relation from the moments in the general case (refs. 1 and 11). However, for 
several reasons, mainly because of the volume of arithmetic involved, we have 
chosen a more specialized approach. We use the quotient-difference algorithm 
(ref. 12) to obtain the three-term recurrence relation parameters from the 
moments. This method will always work when the interval of integration does 
not have the origin as an interior point. The procedure is briefly as 
follows: Define 


. 0 


Generate the rest of the table 


,(j) 


M 


( 0 ) 


( 1 ) 


( 2 ) 


e l 


( 0 ) 


( 1 ) 


^ 0) 


j + 1 

Mn 


j = 0(1)2N - 2 


.( 0 ) 

'N-l 


,(D 

"N-l 


( 20 ) 



(2N-2) 


(2N-3) 

e l 


q 


(2N-4) 

2 


7 


using the rhombus rules 


; (n) = e (n+U + [ q] J n+l) - q^ n) ] , n = 0(1)2(N - k) - 1 


( 21 ) 


(n+l) (n+l) 


q k 


(n) _ «k *k 


+1 


.00 


n = 0 (1) 2 (N - k) 


with k = 1 (1)N - 1. The parameters for relation (15) are then expressed by 


b n = q< 0 > + 4-1 > n = 1 (1)N 
- ,(o).Co) 


with g l - 0, 


§n = q n-l e n-l • 


n = 2 (1)N 


( 22 ) 


Calculation of the Nodes and Weight Coefficients 


The three-term recurrence relation (15) for the first N orthogonal 
polynomials can be written (ref. 13) as the maxtrix equation 


where 


t $(t) = A <£>(t) + B * N (t) 


(23) 


4>o Ct) 


0 

*i(t) 


0 

<f>2 (b) 


0 

• 

B = 

• 

^N-2 (b) 


0 

‘t’N- 1 (b) 


1 


(24) 


and 
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A = 


bi 

S2 

0 


1 

b 2 

g3 


0 

1 

b 3 


0 

0 

1 




b N-l 


1 

b N 


(25) 


From equation (23) it is immediately seen that the roots of ^(t), which are 
the nodes of the quadrature rule, are also the eigenvalues of the matrix A. 
Applying a diagonal similarity transformation to the matrix (25) yields the 


symmetric Jacobi matrix 






" b i 

^2 


“ 




^2 

b 2 

^ ( ) 





0 


- : . 




J N = 

/ 


; 



(26) 




J • b N- 1 





\ 

- 



b N _ 



which has the same eigenvalues as A. In reference 13, 
that if Qi,j is the first component of the orthonormal 
associated with the eigenvalue tj^, then 

it is further shown 
eigenvector of J ^ 




WjN = M 0 



(27) 

Hence, both the nodes and the weight coefficients of the quadrature rule (10) 
can be calculated from the eigensystem of the symmetric tridiagonal matrix 


(26). The version of the very fast numerically stable QR algorithm (ref. 
14) given in reference 12 to accomplish this computation was used. 


9 



Calculation of the Error Coefficient 


The coefficient of the error term (2) can be obtained from the g n . From 
equations (16) and equation (3) it immediately follows that 

N 

k N = M 0 n g i+1 (28) 

i = l 

These constants are useful in obtaining error bounds if one has an estimate of 
max | f C 2 ^) (x) | . 

0<t<1 


Details of the Calculations 

From the formulation of the calculation of - quadrature rule (1) from the 
three-term recurrence relation parameters (16) as the determination of the 
eigensystem of the symmetric tridiagonal matrix (26) with no multiple eigen- 
values, it can be shown (ref. 15) that once the three-term recurrence relation 
is known, no error is incurred except for round-off error. Therefore, all of 
the ill-conditioning of the calculation of the quadrature rule from the 
moments must occur in the determination of the recurrence relation. With this 
fact in mind along with the aim of obtaining a 24-point quadrature rule with 
nodes and weights accurate to 25 significant figures, the following 
calculations were carried out: 

(a) The moments M n , n = 0(1)51, were calculated from the recurrence 
relation (19) in 70-digit floating-point arithmetic. The first a moments 
were calculated using appropriate values of the gamma function accurate to 65 
significant figures. These values were obtained by the method given in 
reference 16. 

(b) The coefficients b n , g n , n = 1(1)26, were obtained from the moments 
via the quotient-difference algorithm. This computation was also done in 
70-digit floating-point arithmetic. Using Gautschi's estimate of the condi- 
tion number (which is remarkably accurate) we find that b 2 $, g 25 had 

k = 65 - 1.53(25) + 2 log 10 (200) ^ 30 


significant figures. 

(c) The eigenvalues and the first components of the orthonormal 
eigenvectors of the matrix J n (eq. (26)) were calculated, using a version of 
the algorithm Gaussquadrule given in reference 12, for N = 2(2)8(4)16, 24. 
These computations were done using 30-digit floating-point arithmetic. 

In order to check these quadrature rules, the identities 
N N N 

D»3N- M °- t29: 

j=l j=l j=l 
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were used. It is our experience that, for quadrature rules on a finite 
interval, these are sufficient to determine the computed accuracy of the nodes 
and weight coefficients. In no instance was the disagreement more than a 
relative error of 10 -28 . 


RESULTS 


In tables 1 through 12, the nodes and weight coefficients of the N-point 
quadrature rules for N = 2(2)8(4)16,24 are given for the 12 previously 
enumerated weight functions. Each of the nodes and weight coefficients is 
rounded to 25S. In table 13 the error constants k^j, for the same values of 
N, are given to 5 significant figures. 

Tables 14 through 25 are tabulations of Mq and of &n> n ” 1(1)26 

of the recurrence relation for the monic orthogonal polynomials associated 
with the enumerated weight function w(t). These tabulations are also rounded 
to 25 significant figures. 

Appendix A is a listing of a double precision version of the algorithm 
used to obtain the quadrature rules. This program can be used to determine 
practical accuracy quadrature rules from the three-term recurrence relation 
parameters. In particular, it can be employed to find quadrature formulas not 
tabulated in this report. 

The parameters associated with each table are summarized by: 


Table 

a 

6 

Y 

1, 14 

3 

-1/2 

0 

2, 15 

3 

1/2 

0 

3, 16 

4 

-1/2 

0 

4, 17 

4 

1/2 

0 

5, 18 

6 

-1/2 

0 

6, 19 

6 

1/2 

0 

7, 20 

8 

-1/2 

0 

8, 21 

8 

1/2 

0 

9, 22 

3 

-1/2 

-1/2 

10, 23 

3 

1/2 

1/2 

11, 24 

2 

-3/4 

0 

12, 25 

2 

-2/3 

i 

0 


In the tables, a ± b denotes a • 10 ±b ; for instance, 0.17 - 6 means 
0.17 • 10" 6 . 


Ames Research Center 

National Aeronautics and Space Administration 

Moffett Field, Calif. 94035, January 19, 1970 
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APPENDIX A 


bo 


QR ALGORITHM PROGRAM OF GOLUB AND WELSCH 


C THE ROUTINE IS ACCESSED BY THE FORTRAN IV STATEMENT 

C 

C CALL GAUSSQ ( M. B* G, ZM. T, W ) 

C 

C WHERE 

C 

C M IS THE NUMBER OF NODES IN THE QUADRATURE RULE TO BE CALCULATED 

C 

C B,G ARE DOUBLE PRECISION ARRAYS OF LENGTH AT LEAST M WHICH 

C CONTAIN THE PARAMETERS OF THE THREE TERM RECURRENCE 

C RELATION FOR THE FIRST M MONIC ORTHOGONAL POLYNOMIALS 

C 

C P (X) = (X - B )P (X) - G P (X) ♦ j= 1(1) M. 

C J J J-l J J-2 

C 

C ZM IS THE ZERO MOMENT OF THE WEIGHT FUNCTION (IN DOUBLE PRECISION). 

C 

C ALL OF THE ABOVE PARAMETERS ARE INPUTS. THE OUTPUTS ARE THE DOUBLE 

C PRECISION ARRAYS OF LENGHT M 

C 

C To W WHICH CONTAIN THE NODES AND WEIGHTS OF THE DESIRED 

C M-POINT GAUSS QUADRATURE RULE. 

C 

C THE CONTENTS OF B AND G ARE DESTROYED BY THE ROUTINE. 

C 

c 

SUBROUTINE GAUSSQ ( M. B* G, ZM , T» W ) 

C 

DOUBLE PRECISION B(l), G(l), T(l)« W(l) ? ZM 
C 



n n r> nnnn nn 


EIGMAX* 


INTERNAL VARIABLES 

DOUBLE PRECISION AA t AJ, BR 9 BZERO» B2» CJ, CT » DEPS* DET» 

* F, LAM* LAM1 » LAM2 9 NORM, Q, RHO, RR, ST, WJ 

DEPS IS THE LARGEST POSITIVE NUMBER SUCH THAT 1+DEPS = 1 
TO MACHINE DOUBLE PRECISION 

DEPS = 0 • 5D0 ** 5A 
MN = M - 1 

INITIALIZE W, SYMMETRIZE MATRIX, AND DETERMINE MAXIMUM ROW 
NORM 

DO 100 J = 1, MN 

G ( J ) = DSQRT ( G ( J+ 1 ) ) 

W ( J+l ) = 0.0D0 
100 CONTINUE 

W ( 1 ) = 1.0D0 

NORM = DABS ( B ( 1 ) ) + G ( 1 ) 

DO 200 J = 2-, MN 

NORM = DMAXKG(J) + DABS(B(J)) + G(J-l), NORM ) 

200 CONTINUE 

NORM ■ DMAX1 ( G( MN ) + DABS ( B ( M ) ) , NORM ) 

DEPS = NORM*DEPS 
LAM = NORM 
LAM1 = NORM 
LAM2 = NORM 
RHO = NORM 
MN = M 

BZERO = 0 « 0D0 


w 


C 



non non 


LOOK FOR CONVERGENCE OF LOWER DIAGONAL ELEMENT 

101 IF ( MN .EQ. 0 ) RETURN 
Ml * MN - 1 
K = Ml 
I * K 

IF ( MN "•EQ. 1 ) GO TO 110 
IF ( DABS (G(MD) .GT. DEPS ) GO TO 120 
110 T ( MN ) * B (MN ) 

W ( MN ) = ZM*W ( MN ) **2 
RHO = DMINKLAM1# LAM2 ) 

MN = Ml 
GO TO 101 

SMALL OFF-DIAGONAL ELEMENT MEANS MATRIX CAN BE SPLIT 

120 K * I 
I = 1-1 

IF ( I .EQ. 0) GO TO 130 

IF ( DABS ( G ( I » > .GT. DEPS ) GO TO 120 

FIND EIGENVALUES OF LOWER 2X2 MATRIX AND SELECT ACCELERATING 
SHIFT 

130 B2 = G ( Ml ) **2 

DET = DSQRT ( ( B{ MN ) - B(M1>)**2 + 4.0D0*B2 ) 

AA = B ( MN ) + B ( Ml ) 

LAM2 = AA + DS I GN ( DET ♦ AA ) 

LAM1 * ( B ( MN ) *B ( Ml ) - B2 ) / LAM2 
EIGMAX = DMAX1 ( LAM1 »LAM2 ) 

IF ( DABS ( E IGMAX-RHO ) .GT. 0 • 1 25D0*DABS < E I GMAX ) ) GO TO 135 
LAM *= EIGMAX 
135 RHO = EIGMAX 


C TRANSFORM BLOCK FROM K TO M 

C 

CJ = 6 ( K ) 

IF ( K .EQ. 1 ) BZERO = B ( K ) - LAM 

IF { K .GT. 1 ) G(K-l) = B ( K ) - LAM 

C 

DO 140 J = K, Ml 
JM = J - 1 

IF ( JM .GT. 0 ) BR = G ( JM ) 

IF ( JM .EQ. 0 ) BR = BZERO 
RR = DSQRT { BR**2 + CJ**2 ) 

ST = CJ / RR 

CT = BR / RR 

AJ = B ( J ) 

BR = RR 

CJ = G ( J+l ) * ST 
G ( J+l ) = -CT * G ( J+ 1 ) 

F = CT * AJ + ST * G ( J ) 

Q = G ( J ) * CT + ST * B ( J+l ) 

B ( J ) = F * CT + Q * ST 
G ( J ) = F * ST - Q * CT 
WJ = W(J) 

B ( J+l ) = AJ + B { J + l ) - B { J ) 

W( J) = WJ * CT + W( J+l ) * ST 

W ( J+ 1 ) = WJ * ST - W(J + 1) * CT 

C 

IF ( JM .GT. 0 ) GUM) = BR 

IF ( JM .EQ. 0 ) BZERO = BR 

140 CONTINUE 
C 

IF ( K .EQ. 1 ) BZERO = 0.0D0 
IF ( K .GT. 1 ) G ( K-l ) = O.ODO 
C 

GO TO 101 
END 


tn 
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TABLE 1.- GAUSSIAN QUADRATURE RULE FOR W(T) = l//l - T 5 


tt =2 


Wfl ) 


T (I) 




w ( I ) 


. 25 56 5026201 7 86422 $46902 3 9+0 
.6712343337313633382217330+0 

KJ=4 

. 7337729738633181679464650-1 
.33017768181 30437138234797+0 
. 7103740067251409197101383+0 
.9642130008208379539527809+0 

N=6 

. 3 504457048024343418449511-1 
. 1 760429136520695590901 530+0 
. 396451 52902 7293972 732 50SG+Q 
. 6449294 58944 7 569740660749+0 
. 8595O5916216510OS24249422+O 
. 98363 7012182 58463 74622039+0 


. 56524283O441481946O830784+O 
. 8369392748839723 1 5091 9407+0 


. 1840204095652859416343037+0 
.3552760729357526979745190+0 
. 4240168147QG0615316641Q73+O 
. 4388668G8124334Q899G2Q09O+0 


. 889331 40 5CQ 9483979 5691 733- 1 
. 1883153151(23423552937622+0 
.252302787 599054859831 4225+0 
. 2829312063 589e5982Q4 54765+0 
. 2936029070014966020373740+0 
. 29909S748759686064Q1GQ6S6+O 


. 5375072026292320269517046-2 
.2810915250282973437079276-1 
.6815424584800657337585396-1 
. 1240874584649518749900841+0 
. 1939Q98213440S96664295713+Q 
. 27 51 1 7097930818073 5066895+0 
.364776003 2207 548QS143Q290+O 
. 4596146246281561328730855+0 
. 5 561 3076407 532476499336C2+Q 
. 65071 80074834021 521628846+0 
. 7390O45472525497O4O93O133+O 
.8199967771198717626798682+0 
n 8882189486S91781674602866+O 
.9418410249417674523104717+0 
. 9787882928700649047395 520+0 
. 9976277673387387424676232+0 


. 13 76992464 7 4Q8Q639Q00Q854-1 
.31 57 508628 72064 46 71 7 7941 4- X 
.4828438326164452830727601-1 
. 6332 9924 199 50048899 2 92826-1 
.7624917075842071955156497-1 
. 8686263 5 56S2 781 54791 51271 -1 
. 951 5522837661942764637529-1 
. 1013183705051350161 506304+0 
. 1056708850057444537581344+0 
. 1085844325240497328488026+0 

.1104228844826455937182769+0 

. 11150514090121461 53691636+0 
. 1120091891610919185913032+0 
. 1 1237O24756310B9O427971 72+0 
.1124858459068593553511341+0 
. 1 125227493785328406014949+0 

tt=24 


N *=8 

.2041981 505827257341810998-1 
. 10461363005751 56032957783+0 
.2443710169441164617315762+0 
. 4210426148229107978965786+0 
. 6102239466795308909477S68+0 
.7846479140515148192768209+0 
. 9182694706544585845538558+0 
. 9906767463370753336401726+0 


.9394755468827131857553192-2 
. 4885880340937446348471456-1 
. 1 172862757957247918O1S4O0+O 
.2104545570059180564392585+0 
.3225648853337444699080334+0 
. 4465428866931218281307671+0 
. 5744120118107212077066824+0 
. 6977755779244984314317850+0 
. 8O0387751632O43433138Q959+G 
.8987540653327876978445911+0 
. 962693595961 571 5229970911+0 
. 9958064 549069209421 871208+0 


. 52O596989337O6115450225O7-1 
.1145639325324057106882327+0 
. 1630713305653662725772673+0 
. 1947610974432407442518330+0 
.2120368987604358993354671+0 
.2190079554326709593010943+0 
.222584661 5448836726138457+0 
. 2232965301 127648869570040+0 

N=12 

.2403632758736658852372809-1 
. 5450676200573899650345546-1 
. 81 70840202372461734435798-1 
. 1042096714922456939030477+0 
. 12136656031 71^9419728460+0 
. 1333825536396215656186495+0 
. 1410957316373476572953743+O 
. 1456124880780362909417237+0 
. 1479949542872908243847265+0 
. 1490970902386087542590827+0 
. 14952335475O452O40O919263+O 
. 1496481092672582819391009+0 


. 2429278828475440309454799-2 
. 1275612399629951258068546-1 

. 31 1 5782404655719027824890-1 
. 5733823454055490258513748-1 
. 9G87Q9S19484Q2929S3240761-1 
. 1312095640297002983967439+0 
. 1776363 19S69977 191 6483667+0 
.2295720646487391432277570+0 
. 28598725046141 1 5078029894+0 
. 3460143313228104664549461 +0 
. 408661 5998702031901413841+0 
.4728885068696724060515173+0 
. 5376223899143700972199642+0 
. 601 77634641 76344056467814+0 
. 6642678888267722075283289+0 
. 72403793239961 124698606S2+0 
. 78QD6964764O68439221O0342+Q 
. 83 1 4067247073703020228529+0 
.8771706348357654398170683+0 
. 916576521 1377072282786839+0 
.9489474005895053985918391+0 
.9737264074306179008332503+0 
. 3904868547992279873 1 77024+0 
. 9989399366660209810408788+0 


. 6229331762389439765950109-2 
. 144019524241818421 5649QeO-l 
.2235205924450900659819572-1 
. 2994073872269130412501344-1 
.3705022319528801995085645-1 
. 4358CE5498883 5 1 1 61 61005899- 1 
. 494 54638133 782575712744S 5—1 
. 5462546105402723985799879-1 
. 5907656164242609514224046-1 
.62822O0S8883855279711O475-1 
. 6590236539227909121104866-1 
. 68377387395384668720O4958-1 
. 7031945723978392905047393-1 
. 7180619634586139568817985-1 
. 7291478556087100795467687-1 
. 7371779959267327551023791-1 
. 7428060234158851816875289-1 
. 7466004105766723379620775-1 
. 7490408914754268788818152-1 
. 7505209773353856705695703-1 

. 751353752030220490261 7649-1 
. 75177886554047290641 10606-1 
. 7519693055044354826765452-1 
. 7520370465441798528666057-1 



TABLE 2.- GAUSSIAN QUADRATURE RULE FOR W(T) = A - T 3 


Ten 


U=2 


wen 


N=16 


W(l) 


. 1M51 7915443 5048844396394*0 
. 7299837021 6O40491O5558436+O 


,662001 66307847 49569696306** 1 
.3137372655764858605014977*0 
.6373331662396621857745465*0 
.9000244758525026819403301+0 


.3265016120507537780640628-1 
. 1636297885906563589690199*0 
.3671361940733304348636674*0 
. 5973299786650457034259124+0 
. 605229969485021 4662385373+0 
.9488328821876800658801781*0 


. 4567144640653860539471959*0 
.3845947991296865027578155*0 

N=4 

. 1656934398410666831255331*0 
.3042133775864010686647449+0 
.2706455105812883441206441*0 
. 1C307 5693 5 1864964607940994+0 

N=6 

. 0281 429551661 712180297748-1 
. 1736068736400336354135599*0 
.2195129152302947223176756*0 
.2012759089066132717405893*0 
. 1256668806154776117658401*0 
. 3843238928623619366436911-1 


. 5227959788538227454999262-2 
.2733674290961945408013557-1 
. 662661 4245331 592251922190-1 
. 1208046761 76286540701 5243*0 
.1883705438556211470321585*0 
.2670962014513314959182705*0 
.3539414243416655439082627+0 
.4457901255296638736811900*0 
. 5393912277383994409714323*0 
. 631 466507 1 7 5623 1 1 42 51 0569*0 
. 7188237864420323069501037*0 
. 7984601932870751 51 5752630+0 
.8676567213358183517432077+0 
. 9e4062903205579e32791 7945*0 
. 9657690685874632Q28O3722O+O 
. 9913698122829605273071696*0 


. 1339271 583690751 303400*3-1 
.3070222991546851968837281-1 
. 46912621 1 5799225324525275-1 
.6137203812169870621885229-1 
. 7344644623057134954268726-1 
.*4884*® 59*5 1187661479-1 
.8786501279812335100577276-1 
. 8904561872582474309630673-1 
.8574122265102774163103146-1 
.7805283886633371764297321-1 
. 665791 578881 1 13859*8831 4- 1 
. 52439896774075242120*705-1 
.3719407307099849880*1719-1 
.2265907115084945477547454-1 
. 1066029973847135256221360-1 
.275752*4282*13552322657-2 

f«4 


Nt=8 


. 1934656167936207757272972-1 
. 99018729223251 72116918954-1 
. 23084754046305639773301 89*0 
. 396922071 7426023399055997*0 
. 5752674877 528963CE 5006790*0 
.7431*1970600495053513068*0 
. 8799190591321317956729686+0 
.96906*76053395*50658543+0 


. 4931331*4303393679641235-1 
. 10817396764031330*758345*0 
. 1514504502617961717881725*0 
. 1704456012399666570353590*0 
. 1591003848655220277664058+0 
. 1193867060920149889416669*0 
. 65066588509*775093968879-1 
.1837224634309571415*7157-1 


h *=12 


. 9056647579583088582247288-2 
. 47t387 5096703297 50*1 68546-1 
. 1 1296864289560*166100345*0 
.202*46303277493996637*1*0 
.3100060545690904280931601+0 
. 4289235601712044714730553*0 
. 5517319234*3368407303446*0 
. 671 104067375860*80142006*0 
. 77996124607954785099424**0 
.87193332*12332*48961670+0 
. 9416890*72880788336906**0 
. 9852104725820194456779*3+0 


. 23 169SQ84 5081 3790969993 74-1 
. *50538335013947873397337-1 
. 784954631 5493191270006718-1 
.3919309759061723038438132-1 
. 1126672241820038586382610*0 
.117149051791496000241*13+0 
. 1115679434129505966942145+0 
. 96267529577355*10502*41-1 
. 7350748512900176678957520-1 
. 47406714962963631 1*75478-1 
. 232246671 1095466042813672-1 
. 61 54794482043808870969930-2 


. 2384317*1299CBQ051288686-2 
. 1251962220725596642954834-1 
. 3057820*396*46018614188-1 
.5626592911546632136424419-1 
.8915875232221724633622034-1 
. 1287129684673059356510456+0 
. 1742749124369632060974863*0 
. 2250927045659366332346948+0 
. 280329928961 669344667881 0*0 
.33908081267999025*188817+0 
. 400386372720338401 1 501618+0 
. 463251 02371 68698591039599*0 
. 5266592377341232458088721*0 
. 5895919693795355817916679+0 
. 6510426573166384947335978+0 
. 710032670*17*111709*11+0 
. 76562608920595190841177**0 
. 81694369670718*106896303+0 
. 863 1 760806964236057497395+0 
.9035955475035686094451520*0 
. 937567123565496489137*51+0 
.9645577917305635933491144*0 
. 9841445229716109399633630+0 
. 996020561 030613*19204399+0 


.6113993099644*1960953168-2 
. 1413431694865300847045019-1 
.2193310967954520235966209-1 
.2936859537300242183907772-1 
.363107**170591508543291-1 
.4263011007850173362339*1-1 
. 4819465411682937572833047-1 
. 5286913461665508335494163-1 
. 565169955903475O141296166-1 
. 59OO667658236846051O12453-1 
. 6022223810226121 1 71817134-1 
. 6007737654356975829878755-1 
. 5853092795436701 176791383-1 
. 55601443861 7435491 7 542 1*-1 
. 5137839027009653252979793— 1 
.460279110202844787*10266-1 
. 39791 77969900672738902021 -1 
.3297899*3722074418537579-1 
.2595037140575991551149700-1 
. 190971 5O9Q08197 57391 49705-1 
. 1281 5547781 12044601631 509-1 
.7479294539053714968893715-2 
. 34127185947974781 1051 1 767-2 
. 8667444736337065125761007-3 



M 

O 


TABLE 3 . - GAUSSIAN 


,2®t669i5S3G3C»Si9ai290«3S+0 
. 06i9?S2674539eS2l3S0Sa23«+O 

>*=4 

. 7277^0*055007^ ^21 522908-1 
. *4eaT3S4S5*442e9t2 5272824+0 
.70611 44^ 1CU344Q1 3278365+0 
. 9«3^133Se3CXH4iQ740e77+G 

f*=6 

. »46 51 73 12 3 3 37 i <2701281 680- 1 
. 1 74938238631 591 363363192^+0 
« 3939074C6»55073«56706476SH0 
. 641 73T9976249533SSS22650S+0 
. 8575«341641243TSe9Q87216S+0 
. SG35630566E21015906D5S488+0 

r^e 

• 2O534961722577847405833I3?*- 1 
. 1041 51 403966721 94908282 52+0 
.2431927120162781 543395881+0 
. 41 SQG841 42307212670043 320+0 
. 6Q77S93&G 503663491 8659123+0 
. 762SS28C23739933Q031 77365+0 
. 9173344734908627201617774+0 
. 99O5582726e4931498S09C323+O 

N^12 

.9368466185410966624422173-2 
. 4871 850369540044820289467-1 
. 1 16931 79532938<23531CKE40>0 
. 2097785777 560890452133455+0 
.3214931648251073067438977+0 
. 44510571 76953937910104021+0 
. 5727818195467991104532969+0 
. 69623331868631 4548301 4736>Q 
. 8071 63 581 9181629481032872+0 
. 89S0O2163241 65957 5 1388542+0 
. 9623889657870023748025564+0 
.3957706104214381936614374+0 


W<!) 


551659S59O4O709Q4O1463885+O 
7 590691181052 703650860313+0 


i {323^17961 1782 363 7 3 53400*0 
343828505 7 5052888902^162+0 
, 3939416^10848728754550179+0 
, 30737O4CDS98875 5 57001 7437+0 


884301307575301372904891 5-1 
18S8iCB38GB8548508313222O+O 
.2459021035337356215535578+0 
.2659593732418354633443651+0 
. 2645479871101 195446S06652+0 
.2595732426742706300401405+0 


. 5184036125414334056201225-1 
. 1139537670579674370472625+0 
. 1613009010587979044010186+0 
. 1895547135053954583285030+0 
. 2QOSOS1 307333642056987556+0 
. 201031 548259448CO93491039+O 
. 197691 5046560574000359^1+0 
. 1 9504493061 08860592397948+0 


. 235586912723820850127 5233- 1 
. 5434176351833615658554056-1 
. 8138173830717595434183212-1 
. 10045320452 38061422463910+0 
.1195493063158403266422300+0 
. 1235722726671372309621648+0 
. 134394S61O105444212658845+O 
. 135519O7158125916653O0667+O 
. 1345486462694183752769665+0 
. 1328069309951351116629661+0 
. 13120669921 50097255780927+0 
. 1302854866679643492262125+0 


RULE FOR W(T) = 1//1^ T 4 


ret) win 

U=1S 


. 5333721599382731540284820-2 
. 28CN89692444579031O399a3 5- 1 
. 66aM!42765Sea8279©779930- 1 
. 1238024685433812891081035+0 
- 1S3443331 952807872 56282 500 
.2744319170168170091744477+0 
. 36386273 1 406051 0373066440+0 
. 4585077287258034561 530680+0 
. 5549113307976360899691061*0 
. 6494932643939194555662962+0 
. 7337O5543190674939553O832+O 
. 8191141454896716610300751+0 
. 88760536787997S28361SS21CWO 
. 941494895o507P87S362323S7+Q 
. 97^552953336219753880589*0 
. 997612508S234S4G4S81 50586+0 


. 2425855359029000582910457-2 
. 1273793657151042334049450-1 
.311128S467700ED2621255118-1 
. 5725401491664381 K761 5139-1 
. 90734122 549S81 32762967 574-1 
. 1310059649837677896060252+0 
. 1774118522O1120DQ53113855+O 
.2291942266157932604209415+0 
.2855077134260338168776594+0 
.3454316362030141706420780+0 
. 407993301 9S332Q5772 17 1244+0 
. 47213203571 56853965573961+0 
. 5368142406568784689392195+0 
. 6009496943254136107066142+0 
.6634590661594063194261164+0 
. 723282 53622405 585001 6 539Q+G 
.7793986155066144429755630+0 
. 8300432729683367373655738+0 
. 8767280636397445933818945+0 
. 9162571307971588450936636+0 
. 9487424866297988725038430+0 
.9736172284448392699205196+0 
. 9904464126102434692363390+0 
. 998935379003045757 5600607+0 


. 1 37407 588860407763 £»9744 7-1 
. 31 •05Sm85729C66970a46eei -1 
. 4616355439644947104238779-1 
. 6310431399806402976285 S22-i 
. 7582072337587391 1CE221 563-1 
.8596799165278574321209731-1 
. 934206100961838971 5690461 -1 
. 9031 lC©12677a3S2CE77Q5444- 1 
. 10inaSO^?7Z37 476364603831 +Q 
. 1 0200429976330594 56859069M3 
. 1018591438194 1080226 55226+0 
. 1010626303693790566053031 -KS 
. 10CD121G2S37G9846339324S2+O 
. 99GCE4 381 7370808274630027-1 
.9823110965381250431786479-1 
.9781757777520341522264673-1 

N=24 

. 622D49Q1 564 54 10S3 5330461 5-2 
. 1438 12493290201 5834345755-1 
.2231 881 C3Q 5663840954 12982-1 
. 29eae2 77484237 5 5431295 501-1 
. 3697799066248393004834328-1 
.4346356351924814633463833-1 
. 4525250727419372968532237-1 
. 542704S31 7787278124923864-1 
. 584723599457311 5090330080-1 
. 6184793122028814251605964-1 
. 644239891877492101 5276560-1 
. 6626187866854466123920923-1 
. 674506S24179354896231356O-1 
. 66O97O9585959281C09CES665-1 
. 683 1 471 121 799E3 1337891 495-1 
.6321469633 76 59ei3 54974691-1 
. 67390622462455879423841 53-1 
. 6745 554813391 492589377063-1 
. 6S95S3S3 544488594077851 53-1 
. 66466055341 53252707928378-1 
. 6602334883927204500440237- 1 
. 6566559244425246597530933-1 
. 6541444586926673242177174-1 
. 6528530366620744722669747—1 



TABLE 4.- GAUSSIAN QUADRATURE RULE FOR W(T) = /l - T 4 


T(l) 


. 1 $> 7 423 794 9 1 6 7 50 5963 1 59 729+G 
. 736326088 548995 551 4029448+0 

.6671843902813787620345253-1 
.31649O87365QG1437251600Q1+0 
. 641 4706667330448965993946+0 
, 901 85082 7 533063 7907392523+0 


. 32824491 45482232 564997399- 1 
.1645854675152141536378311+0 
. 36937952478489CO463213310+O 
. 600281 7 72 1 64864864 743 57S4+0 
. 80741 1 58899053 57861359224+0 
. 94954461084461293329239100 


I ) 


.4655093301964751153818097+0 
. 4085098545675648214398035+0 


. 1670728285270703586380140+0 
. 3101 6492 59e031363363437G9+0 
.2859274812041764659755736+0 
. 1 108539491 124794785736546+0 


. 8326 70555962280 5018S08282- 1 
,175041777 5079680941 773229+0 
.2242074446510954033681721+0 
.2116428118176476084722913+0 
. 1368097432847228295498204+0 
. 4297035190637795106492362-1 


TO) 


. 5238880081916714983811427“2 
.2739459335892662143545982-1 
.6641016141603248718859137-1 
. 1208773707265071 554448609+0 
. 188815572353249326961 5684+0 
. 2677 5034905426391 843841 73+0 
. 354S1C262GB77389126785167+Q 
.4468458635329833140541219+0 
. 5405624097281343646571991+0 
. 6326524723605557412259872+0 
. 71991691462677761 59266629+0 
.7993695451731488025810808+0 
. 8683262 50234027864161 6361 +0 
. 9244809382923489334656035+0 
. S659S95 553623580566310523+0 
.9914218919211441955502528+0 


u i! 1 

T4= is 

. 1342077082401426204394644-1 
. 307686320584948 5 1 663 7 1 82 1 - 1 
. 470e 5832836 5401 5 1 9e3 7 7642 - 1 
. 61 57021407626646221 132632-1 
. 738386330876613862441 7644-1 
. 832885 128683 5 1 42 7 7 5829836- 1 
. 893 72 5623 562 7862 5 58 583 7 59- 1 
. 91 5491259651001 1397785308-1 
.8937068074492539138415139-1 
. 826561046363399741 43 70908-1 
. 7169556718806446708912261-1 
. 574009136 197056361 0952680-1 
. 4132146848204729345341 804-1 
. 25487935 102273403266624 59- 1 
. 1210334840869214829096594-1 
.3148882509284341675161869-2 

N=24 


KE 8 

.1942520141529435294409544-1 
.9944388251406683185069363-1 
.2319197278934014769457291+0 
. 398771 848289871 7747730782+0 
. 577 5420360065801 57746C006+G 
. 7451 5603682 5 596 524S072 726+0 
.881 0939400075294263991 941+0 
. 9694085376125641 6S23471 53+0 


. 90816252432 5672289SC68787-2 
. 4722051495678480760727549-1 
. 1133C3234 7 146602956248536+0 
. 2031 5939107903233377 59724+0 
. 31 10893866019913446425329+0 
. 430272 54241 8681 98S92 54 738+0 
. 5532808759994113664281759+0 
. 6726227736866722132128896+0 
. 7812206713912676424811427+0 
. 872787666475724377 58S091 4+0 
.9421199469616726476371718+0 
. 9853263381 762 561 1 6891 9370+0 


. 49516249462334O4137O05629-1 
. 1087190937556033325146374+0 
. 1 5294131323032001 57670059+0 
. 1744820864755212398896874+0 
. 1667782489385567568343932+0 
. 1267989453709585213439622+0 
.7208635345043414406078184-1 
.2071 688842971 188504108901-1 

. 2323413935226197981026063-1 
. 520S1G9343O1 5637064046647-1 
. 7879716147785233806827536-1 
. 998S68573582S6CQ275083332-1 
. 1 1422371 04171352888480175+0 
. 1202414402703363646769067+0 
. 1 165272169528736613666480+0 
.1026418157218367640555386+0 
. 800472680132229395263 7477-1 
. 5261557786525529613676088-1 
.2616351966249146671631370-1 
.6999377682250864225211225-2 


.2387672779195479726520233-2 
. 125373393 5563277575134309-1 
.3062196072467974986062373-1 
. 563478S1890975O9967O73675-1 
. 8S891 72931319036671696595-1 
. 1289105406841 53455933091 5+0 
. 1 74 5 50606 126G063 189 584691+0 
. 2254584382055001839489014+0 
.2807937060964253575960866+0 
. 3393442192m 8204771420740+0 
. 4O1O426143953767O0564788O+O 
. 4639841 5621 77834398552379+0 
. 52744A9095439605545751 563+0 
. 5SO39957690S48S6311946681+O 
. 6518384787695364557562333+0 
. 710783395861291962820661 1 +0 
. 7663Cei8941 549ei3775130GD+0 
, 81 7 52221973 1 91 97284984807 +Q 
. 8656426023138325744672855+O 
. 9039453087652278701321 1 10+0 
. 9378050723789461462899838+0 
, 964696019043099953 5300406+0 
. 984208922261 8837 57033448 5+0 
. 99603697^784407389501067+0 


. 61226073-10795912904555501-2 
. 142 544771 5261 713974120604-1 
.2196541888346776646552949-1 
.2941 554037O453583958729O8-1 
. 35380237 542195671 1 7567382- 1 
. 42741 19693389322756040533-1 
.4838476058695132559350565-1 
. 53 1 93261 S2693 54462 5 1 10906- 1 
. 5707187329531747679969179-1 
. 5968690613527000043672351-1 
.6152804689496367183115263-1 
,6188654673820973947939941-1 
. 608759130378861 5205899223-1 
. 5845167940738630638284983-1 
. 54632862 5539G87Q 52861 1835-1 
. 495209476479775231 5424432-1 
. 4331 18981 5247003 5 7 4843679—1 
.3629792369879171083324069-1 
. 288571991 79387 122 5873 1 853- 1 
. 2143l639e 513 5783868440588-1 
. 14494638597819O605413O7O4-1 
. 8512 1 8270802476 58 71 4 1 031 7-2 
. 3901 732853641 61 87 572381 92-2 
. 9936884534980447383697103-3 



ro 

N> 


TABLE 5.- GAUSSIAN QUADRATURE 


T'l) W'l) 

U=2 


.2275715979585846041106961+0 
*6571031 565242673380969259+0 


. 7212564772417261791344756-1 
,3432544710951126068353072+0 
. 6996260160266079844668651+0 
. 961 81333601 7496945371 5734+0 


.3464064059329285167774036-1 
. 1 738371378996293513536924+0 
.3910263221669966170957644+0 
. 637335704127966656341 5256+0 
. 8 544097924096136674234567+0 
. 9828626764456324881865124+0 


. 21E41 42560063658862619342- 1 
. 1036566185863387993093046+0 
.2419436207017047237525148+0 
. 41665734O129QO2987O433703+O 
.6045203635673331562034694+0 
. 77966281 54694230593706342+0 
. 91 57667296761967301394336+0 
. 9903514267781723841627602+0 


. 5396250766377089967178524+0 
.6747002473060618091921185+0 


.1807113803650086784537596+0 
. 34033557343662105699101 57+0 
.3645092051233290614438561+0 
.3287691650188320090213394+0 


. 8780812226009906562400217-1 
. 1 851 97800E91 3 5 599009533 12+0 
.2413127540563494242581581+0 
.2504770241229544842489740+0 
. 2327459354897027858484999+0 
.2167036877233290558350055+0 


. 5160038594556015328914058-1 
. 1133847276574747Q03465534+O 
. 1600835945713090655364412+0 
. 1859291321 5SS3Q3 10601 7261+0 
. 1906568817477191072031545+0 
. 1815639112785202922248961+0 
. 1693780894239286477538975+O 
. 161728801 162648528954161 5+0 


tt=4 


N=6 


N=8 


N=12 


. 9339328521 134557282451370-2 
. 4856504096397518720138830-1 
. 1 16552813783338141 5563990+0 
. 2090631287133629190469690+0 
.3203219461689760718843399+0 
. 4434076804369578485004445+0 
. 5706512113330910886379781+0 
.6939798832691580721760148+0 
.8052557693724250451817117+0 
. 896762781 5061241 706578995+0 
. 9618689689121458455218392+0 
. 9957063921 729528527991001+0 


. 2389393891292663344950964-1 
. 5416747728559574390458385-1 
. 8109567799504302844644799-1 
. 10296051 507511 51068775399+0 
. 1184648909956516431709682+0 
. 126970187789237898851 1 544+0 
. 1288947829218173402447284+0 
. 1258824054659614489728578+0 
. 1203155242449796183664725+0 
. 1144482159786554418766757+0 
.1096515170235062445182999+0 
. 1073801902553006552307126+0 


RULE FOR W(T) = 1//1 - T 6 

T (I ) W I) 


. 5351113809427072606442839-2 
.2798258465146087668230399-1 
. 6704O91218342635634521OO5-1 
.1234979286140205703158452+0 
. 1929491916472123105349557+0 
.2736960922637372392879847+0 
. 3628396297576300897779148+0 
. 4571847529093397375210405+0 
. 5533389146000698023264119+0 
. 6478052277124858827437712+O 
. 737O731O57314719776452906+O 
. 817730197621 5492070672670+0 
. 8866027466593136263792395+0 
.9409124766791088735150425+0 
. 97842729681 546O424762606O+O 
. 9975861 131075394273035640+0 


. 2422000286060772053054837-2 
.1271775477525497549545132-1 
. 3106331059803663470142352-1 
. 571619655334O06367O59368O-1 
. 9058629550459076955082901-1 
. 1307885185320622863338079+0 
. 1771099655810435804983677+0 
. 22879021906478Q58485Q38O9+O 
. 2849839071662064658939877+0 
.3447847755062304653840692+0 
.4072004446635404844282340+0 
. 4712171 585332585916564732+0 
. 53578605263O1514757894338+O 
. 53984299CE290D44412233 525+0 
. 662322601 5787848020540519+0 
.7221734173631736578407328+0 
. 778374703288243567951 5398+0 
. 82995493800374797441 11974+0 
.8760112162319662517110063+0 
. 91 57285004212557257250008+0 
. 9483976360616956069622984+0 
. 9734312573741479128495390+0 
. 990376981 7821039656283918+0 
. 998927524132 562488103 5366+0 


N=16 

.1370841062528893609173151-1 
.3143065893946315022316427-1 
. 480463814452359299581 7756-1 
. 62932389581 1 7820393657079-1 
. 75554181O5775377898406364-1 
. 8547464622238238230948059- 1 
. 9240197342630746434482978- 1 
. 9626708794838849642406876- 1 
.9730415879197964684931308-1 
.9606618044443370154183532-1 
.9332540711199260586974823-1 
.8989554381515707508335647-1 
. 8647291336496365034536576-1 
. 8356202379O972O4O491 18719-1 
. 8148061961 589574737035017- 1 
. 8040274778239699608672771-1 

N=24 

. 6210656995012059540275417-2 
. 1435837824878390925263551-1 
.2228288200453923518582349-1 
.2984305144231008558284727-1 
. 369141661 5902560480010051-1 
.4338043472611628637824633-1 
. 491367Q9O86278555590687O5-1 
. 5409130446950003302082055-1 
. 581 70245681340601 57939301-1 
. 6132362996607207633886CE2-1 
. 63534O26834854O02348387OB-1 
.6482516855195177557517142-1 
. 6526747148750469364294950-1 
.6497665266318900586911039-1 
. 6410377976721 584218808820-1 
.6281865117725849228626807-1 
.6129118953724405826466743-1 
. 5967582627610940706470057-1 
. 58101 86850653308137829202- 1 
. 5667018278654222346787269-1 
.5545466539238632231591124-1 
. 5450634773204470107835428-1 
. 5385825842902753968113112-1 
. 5352977602294919347137555-1 



TABLE 6 . - GAUSSIAN QUADRATURE RULE FOR W(T) = A - T 6 


TU) W< I ) 

W=2 

. 20059979091 699934324001 13+0 . 4740991886530853 343040270+0 

. 74 337266967301 2373 60 10423+0 . 4366448043048177701284513+0 

*=4 

.5709005978379390283597651-1 . 168 54470S32 1499722 79 59661+0 

.3195587357275720163337610+0 . 3150404130016933665797106+0 

.6472033373752203309482565+0 . 3022174823903069664484664+0 

.9047467142100877677667062+0 . 124941 38Se 443430486083349+0 

N=6 

.3301885643412491154815001-1 .8376507044377595806290145-1 

. 1656148590788018079626461+0 .1762928849233464878475253+0 

.3719440513073763233104491+0 .2278783213265854585495529+0 

.6O40534252879632473644538+O .2218957826888740290687443+0 

.8107505487088011779818326+0 . 1 51 1573QS98S5 170104103906+0 

. 95070437283367491 52835526+0 . 4975462659074335969289359-1 

N=8 

. 19513QS825743O02864735993-1 . 49741482985053231 55883309-1 

.9990659071681994909911413-1 . lOSe 51 56O6C09633 7 184021 18+0 

.2330751160369460966485501+0 .1540142321790286987612533+0 

.4009198967214111018079766+0 .1774132739655488599103919+0 

.5805491864258133907482330+0 .1739923449628038957699110+0 

. 7481094100782634165877557+0 . 1400542623555054265461363+0 

. 8829642 543 5449631 057823 1 7+0 . 81 S3 509349243441 38481 7405-1 

.9699753665725955758262449+0 . 2434 1 742 197 5 1CK061 97 56674-1 


.9109553248461558923743778-2 . 2330578289775986737083279-1 
.4736745351176233563132086-1 .5282769168009784050521913-1 
.1136651007748527051514341+0 .7306816327864060351218903-1 
. 2038391 4007948821 79634748+0 . 1 CEB2 13020375723242603937+0 
.3121950750777510842478183+0 . 11518191O576O9958SI349422+0 
. 43 18S84 530958 1906 3497 6739+Q . 122440533233CB8249 5733843+0 
.5552889358442620718172042+0 .1208974487476286135419935+0 
.6747709622272935615857870+0 .1094847784057625785392580+0 
. 7031302393106788384406737+0 . 9826270654202248466544926-1 
.8741483751766431951967678+0 . 5996132602638993110073569-1 
.9428081688576118996431721+0 . 30645461 84885796213196887-1 
.9855201220048573167842806+0 . 8346885633964064205591725-2 


Tin win 

N=16 

. 5251O9193QS927924945165G0-2 . 1345209969399729269995129-1 

.2745883091369993856941566-1 .3064142723857861372194006-1 

.6356802214227931184171355-1 .4714062576220562967593744-1 

.1211713538501204050099609+0 . 61 734906953 795474963001 76-1 

.1632912399212662306598531+0 . 7408966568087104073737401-1 

.2634561183947082689525087+0 . 8374129000188018533748346-1 

.35578795716S2728275457O07+O .9008306597057509300556915-1 

.4481079470745539638190390+0 . 9332400812810480739038738-1 

. 540O54744O9S787372S63136O+O .9044881495738857576420798-1 

. 6342828457076905365824 79e+0 .8706047680708909000931296-1 

. 701 313712971 5631261045001+0 . 775835660751 5352393935734-1 

.6007647661971971503361267+0 .6379009737986536878468816-1 

. 6693926365905363634744017+0 . 4713816366S6O18225O578099-1 

.9251648851214914186909425+0 . 2974742325663 529e667064e9~l 

. 9663021 524482254446210369+0 . 14377706S6904S609e821332S-l 

.9915095220901535354093542+0 . 3782410091032051299462033-2 


.2391423526684103488920106-2 .6132251864223953388690883-0 
.1255706973557172598019554-1 .1417685538270603940568587-1 
.3067O462O79O7868337713967-1 .2200055167170374120707319-1 
.5843787845169766420527154-1 .2946362085805960301063262-1 
.8943619453818105740777429-1 . 3644242691921 109327537358-1 
.1291228594367632265857156+0 . 42821 80228837^60522 50396-1 
. 1 7484497 K3927 52 1699347984+0 . 484959839442621 4627468S1 8-1 
. 22 584S83 1647763 1190904219+0 . 5336919450044786808646035-1 
.281296S9079658990731 13144+0 . 5735469310093953585084064-1 
.3402719281384675172097809+0 . 6037155887001235743291006-1 
.4018013633538108967836668+0 . 6233936531848027886093251-1 
.4648704320680175884990508*0 .6317272633679335716748727-1 
. 3234416187607439187605798+0 .6277982647384900529829680-1 
. 591474854883829981 1226788+0 .6107000703801017445096946-1 
. 6 529476431 801 80E21 1372656+0 . 5797403580363990042497086-1 
. 71 18739033161202337489638*0 . 5347580606068904485390416-1 
. 767320401 6S271G0758446176+Q . 4764997009037507205974737-1 
.81842O5453886782733735780+O . 40694539301281 73643267889-1 
. 864385 568799087201 8S026 12+0 .3294847954909462393818544-1 
. 90451372307900784685461 50+0 .2488491545665348761515850-1 
.93819793S889029736SO5S241+O . 17078836691 7066739Q5S5913-1 
.9S4S322830534519S1 8639740*0 .1015052500651794811795453-1 
. 99431 73S5324d35$552901910+0 . 46940082787 5S73 5393849997-2 
.9980647578184413762566587+0 .1201958853848830877098654-0 


M 

CM 



M 

4 ^ 


TABLE 7.- GAUSSIAN QUADRATURE RULE FOR W(T) = 1//1 - T 


Til) W(I) 

N=2 

.22 339*06 7 1093922 5 3 496 523 1+0 . 5339167766881149659655497+0 

.8499161833327369986334796+0 . 6296757925301 544093369684+0 

t*=4 

. 71 76293263900066991529586-1 . 17978811313O3277705561749+O 

. 341332311 58996667922 573 16-M3 . 3377840978239945065155174+0 

. 695236239930197121 7720250*0 . 3512912563439316855275457+0 

.9604691020524162648196490+0 . 2947291039200154127032801+0 

K*=6 

.3452139477390534675785993-1 .8758359092491355523246647-1 

.1732177736598475398829786+0 .1844983465201171375066967+0 

. 38947471 77630029138281097+0 . 23974941848616405641 14986+0 

.6345412420716818955136444+0 . 2442527628940114022470505+0 

. 851 92949843 591 461 02601 633+0 .2159892307850513821211812+0 

.9824606676647490474211754+0 . 1915192216080118417816248+0 

N=8 

.2018838682681026585586613-1 .5146468835733052233122436-1 

. 1033803072806932965040906+0 . 1 130749545770275298833705+0 

. 241 2692 52265340E32449533O+O . 1 59576312685376301 5016761+0 

.4153774454938662682753184+0 .1848037068295133988926514+0 

. 6025035981 1 57713805308142+0 . 1869476691754207981700122+0 

. 7774563396516134566453904+0 . 1721 5640723689182O4384287+O 

.9144718723516904413655563+0 .1535395861210096977005078+0 

.9901704416964351778327957+0 .1420292442356993083846471+0 

14=12 

.9322753745726375675824757-2 . 2385146651163786116292117-1 

. 48478239621 50604459924433-1 .5406976019336749636371705-1 

.1163413973609962524481794+0 . 8094424874703764950963775-1 

.2096731585957883104083876+0 .1027496385618730094368070+0 

.3196931852592804476105785+0 .1181341055436329574150682+0 

.4424702932471916266180862+0 .1262399319742800631109162+0 

.5693716812578721598963450+0 .1270343859061349938011629+0 

.6924735180409074240565064+0 .1218086774116101643310283+0 

.8038305663383317194531057+0 .1132185990026675543411080+0 

.8957609397537849754497026+0 .1043262612247306037582979+0 

. 96 142729CE 52666369213061 9+0 . 9744006342540497916343516-1 

.9956542896153438813700687+0 . 9377543271589204290841865-1 


T ( I ) Wt 1 ) 

Kf=16 

.5343932029438403261691341-2 . 1368999609953713894610340-1 

, 2 79448826 72 4694694 782S290-1 .3138810048071 185206430814-1 

, 6774881286330022940703297-1 . 4 798Qe 5344 79692 539 7 590563-1 

. 1233281490855831276523778+0 .6284307430852471351199931-1 

* 1926784818404740875679212+0 . 7543936210046509286469272-1 

.2732993745305595406772107+0 . 8531818727487739050700088-1 

. 362288G3 747701 78991 594542+0 . 9213948788893768427280534-1 

.4564433911227921186208830+0 . 9572101837156037163811901-1 

.5524064980990168417994245+0 . 9612713549713257604248359-1 

. 646714347365127971 4280540+0 . 937772372481278861 5322739-1 

.7359292685280644784435120+0 . 8947141936531675007382724-1 

.8166860060514728832961848+0 . 8423872675160439632233471-1 

.8858010601406381776639632+0 . 7907048907206236539766873-1 

.9404271215434436894337888+0 . 7471515962031458018250882-1 

.9702322621545105582280394+0 . 7162974293151042832992681-1 

. 9S7 5632481364243 534900224+0 . 7004318055961689501960841-1 

M=24 

.2419817710567806845254686-2 . 62050503758064001671 51 148-2 

.1270625882424325456557334-1 .1434537145726216283879411-1 

.3103514199177617681273783-1 .2226256397969507127467375-1 

. 5710987882620168021 160333-1 . 2931 555613827484820709854-1 

. 90503182031 53569674979254-1 .3687961554553926778181727-1 

. 1306673694391977062350966+0 . 4333879302820568130161740-1 

. 17694373221097OO52278831O+O . 490872896873541 50633CE864-1 

.2285735324844436076723154+0 . 5403128585679183184222834-1 

.2847107002420195408653695+0 . 5809105372007371524465046-1 

.3444360673302032215919102+0 . 6120397945115497146569145-1 

. 4067729834038955206033572+0 . 6332966774736290224867249-1 

.4707033204415411610017633+0 . 6445830312751339387675432-1 

.5351836501402536091894047+0 .6462199433106763119791951-1 

.5991598728838550647189528+0 . 6390571116623033234126172-1 

. 661 580214761 1035757757371 +0 . 6245144426564912471095966-1 

.7214070972128249941041801+0 . 6044971076956346216627875-1 

. 7776296296442802816990275+0 . 5811 796632 19ea2873 5088349-1 

.8292785486413161546346536+0 . 5567258591729055552211006-1 

. 8754443146924332013451420+0 . 5330424284587717574830342-1 

.9152975318937119960598474+0 . 5116363275836906163098435-1 

.9481098532133713975802415+0 . 4935866602664662503902905-1 

.9732734250847837980318468+0 . 4795963406133997612562349-1 

.9903174112644428881127237+0 . 4700916405203598339433257-1 

.9989207487226400543359058+0 . 4652900858721676319265416-1 




TABLE 8.- GAUSSIAN QUADRATURE RULE FOR W(T) 


A - T 8 


T ( I ) 


. M2375O0853«2O43aa2O0O76e+O 
. 7MQ391 9231189O5707124972+O 


.6762020646034779647 513202*1 
.3212534163069693601322000+0 
. 69099106760234031 55216055+0 
.9070043041796462032666609+0 


. 3313027475666341321 163667*1 
. 1661095510552647647412969+O 
.3733520725141671633334912+0 
. 6067400023246180999623601+0 
.0132212906592681773150661+0 
. 951639397 5671376047068676+0 


. 1956344262078435102748306-1 
. 1001683179778967273816619+0 
.2337O93202421O58867791785+O 
. 40ei046109961 527789666053+0 
. 5823926129399936559551620+0 
. 750194493031 5226764644000+0 
. 8644225594638400688269185+0 
. 97O43903846849643694O9855+O 


.9125565122275229627343619-2 
. 4745121654476581940468389-1 
.1138687866611254531722562+0 
.2042136779918064104489568+0 
. 3127954783581012200101 503+0 
. 4327561644352037276574994+0 
. 5564929629573438845021 162+0 
. 6761932480091 1 70977483304+0 
. 7845189000322064234013185+0 
.8752077517796296135575068+0 
. 9434040098971910195901 522+0 
.9856813991375511742739003+0 


W(I) 

. 478554393840353679979O999+O 
. 4523196631342618202629146+0 

N=4 

. 1 6937638234 59224936299348+0 
.3171120744648610639090086+0 
.3099764357293958421252360+0 
. 1344O91644344361005778351+Q 

**=6 

. 8404924 1 912621 5379C036751 -1 
. 1 76934286584121 5075658857+0 
. 2291 773225223383328082682+0 
.226244O79O5O1443410582448+O 
. 1 597483024602636474666288+0 
. 5472082444512613344291946-1 

*4=8 

. 4987029123993482666537866-1 
. 1095438518452178839518468+0 
. 1 544817O4507462Q83O549787+O 
. 1 7836639800623067 59258848+0 
.1768156337802232964489838+0 
. 1461019285986481288694083+0 
. 88586341 143 134OOPO0CI3O479-1 
. 2710790785376460524522864-1 

*4=12 

. 2334679306084943762837985-1 
. 529218935671 5093855202448-1 
. 7921359243826427639806870-1 
. 1005217345566718070102287+0 
. 1 1 54864769656860636827O39+O 
. 1230779583387508235723436+0 
. 1224393795865288369901904+0 
. 1 1261 79475CB5951 04429821 7+0 
. 9295187663587300418967146-1 
. 648885 500233628532 1680844-1 
. 339918663 5382009889602390- 1 
. 941 5987925061955675749333-2 


T ( I ) 


. 5258O0686423312O5592O1 1 53-2 
.2749557016350563490927281-1 
. 666576397143745920031 1490-1 
. 1213363831011425383909180+0 
. 18955391676786061 51 586221 +0 
. 268839977083649541 4400037+0 
. 3563 196451289120420621 755+0 
, 44881400951241 40749371 1 13+0 
. 5429 5 54 5422796 103081 001 00+0 
. 6353254032171 594767530392+0 
. 7226164732498821281489782+0 
. 0017948431 57056943251 7410*0 
. 87022191 56123662880283397+0 
. 9e57169CE3 530737913066463+0 
. 9665772142195596758646960*0 
. 991 582989901 5769094809478+0 


. 2393570728835934495465314-2 
. 1256838266854397801943109-1 
. 3069812789406790296946786-1 
. 56489CE 164668 579476 112148-1 
. 8951777034604105274728961-1 
. 1292417016676006735154758+0 
. 1750079172634822436804181+0 
. 22606395836774593471 71 650+0 
.2815699755360449733519266+0 
. 3406122049456964283938684+0 
.4022176008435241666455611+0 
. 465369679202 5484559462736+0 
. 5290259O4148542O875628999+O 
. 5921370084927656173952666+0 
. 6536681306503592327856844+0 
. 7126205920123813369673510+0 
.7680522014950836521114354+0 
. 8190940D53883210571 707650+0 
. 8649624961387718617281908+0 
. 9049677505743729235272227+0 
. 9085186698O6O53O35O547735+O 
. 9651266397406279838707362+0 
.9644063668641725384371393+0 
. 99608821 G3CE86633 1961 6230+0 


wm 

N=16 

.1347003381 579255481 527408-1 
.3088284802541602794364957-1 
. 47205098265449244 50457275-1 
.6182151044238530190608770-1 
. 7420051616977847969688565-1 
. 8388976842622197297529103-1 
. 9O52494278522999726970196-1 
.9381115409519722873095513-1 
.9346690447570131313307361-1 
.8915656592941315824227053-1 
.8053057320649268033330002-1 
. 67571 558750166297327 36073-1 
. 5102490594583470374182135-1 
. 329660986CE73O777957ae424- 1 
. 1614968672176494128266399-1 
. 4294O96317O4O02O58CD1418O-2 

N=Z4 

. 61377399O075O2O45O523O539-2 
. 1418963084148020463342254-1 
.22Q2O50O3572O592O78733463-1 
.2949060120632673188358820-1 
. 3647629974569800405851335- 1 
. 4296256673848763307263197-1 
.48544228272421 85 549232 51 2-1 
. 5342745042256629671156613-1 
. 5743O729402416829O749656O-1 
. 6048467770287235681060999- 1 
.6252950454379717638162151-1 
.635091774774436137916561 5-1 
. 63362805251 5336001 5031298-1 
. 6201 7250563258871338051 73-1 
. 593889401391386330571 5867-1 
. 5540577482982729789077672-1 
. 5005878750630705283204224-1 
.4337856431691832145559328-1 
. 356689GD77889422388875 182-1 
.2734596744591389112474242-1 
. 19023944941 56246980952381-1 
. 1 143 5541 80223441 764093889- 1 
. 5333525376596788243009912-2 
. 13729469029611O5627644295-2 



N> 

G\ 


TABLE 9.- GAUSSIAN 


. 13030212209544 72145906203+0 

. «3529407216t38547308672S9+a 

N=4 

. 5 367 123 592*3210283287666*- 1 
. 2939560710570540734322316+0 
. 676301 780650064268481 1674+0 
. 9 593 732452 797412837682852 -HD 

^6 

. 1630137480570383040237659-1 
. UQ78971E61 48294074386980+0 
. 3 59746901 8S26891 332442995+0 
.6183077970552631336327448+0 
.8477627723245529340650405+0 
. 962193970920439905201 6066+0 

NJ=3 

. 9291941626561905519963069-2 
. 9166048171091 533857620968-1 
. 2162471487736233373723057+0 
.3940389663179270962833099+0 
. 5891269784792185760052748+0 
. 7717836561446507583624536+0 
. 91303331 54316830469500974+0 
. 9900678400662 4089944334 59+0 

N=12 

.4192797654613511262829163—2 
. 3723863 509774276958814549-1 
. 10121 56291 Q6953CB71 740544+0 
. 1919896451304628650909087+0 
. 30368765781 55392305757587+0 
.4289917621729014214568807+0 
. 5595311070683353676119461+0 
. 6863967833O9571O063294938+O 
.8007653032146166233036092+0 
. 894565727553770751 5962959+0 
.96111 102184681 1 5377298035+0 
.9956263595197316077650052+0 


WU5 


13Sei29e83 523042851 524713+1 
1036521364364538760295228+1 


7471880363857365358553503+0 
6588793946745730506333283+0 
. 5431246934004264188603391+0 
. 4776585234268436062703239+0 


. 5O0O762231947474424783142+O 
. 4777325646804528407335396+0 
. 426281 7602419202144594560+0 
.3725351435282757707329870+0 
. 3322997355697048849407952+0 
.3117232206724904584690497+0 


.3844316310299928436073746+0 
.3709115736153533562649943+0 
.3458423218533769927317627+0 
.3143703092123741543738179+0 
.28328559SQ920694622058603+O 
. 257795865521 7896138176260+0 
. 240367090252476251 7160444+0 
.2316462553101489371024615+0 


.258347163209365857331433O+O 
.2541602290011967910205067+0 
.2459837928299936620355372+0 
.2344101057468472728018143+0 
. 2205448139203199875780012+0 
. 2O58278948127O05QE6687837+O 
. 1916574391246276673771 197+Q 
.1790897087931024703817512+0 
. 1687569506336595171441316+0 
. 1609502704813483415365915+0 
. 1 55755O0634O9237765498344+O 
.1531671929634917653944373+0 


Til) 


r*=i6 


Wi!) 


.2367407232573675923254878-2 . 1944733719926803951039243+0 

.2117526341205507221766248-1 .1926752464149563768611977+0 

. 5809629413084753512523175-1 . 1891104744182135918024864+0 

. 1117690621256004629234749+0 .1838979877217325404372531+0 — 

. 18G21 5648G34558197651 1365+0 .1772677547308297307861596+0 

.2600075346745081337067323+0 .1695840313372804638882497+0 

. 3508443714049945296521878+0 . 1613066250140909026470397+0 

. 4466481420507367879554784+0 . 1 K91 8628354 597098S2S6469+O 

. 5446767586204551688184804+0 . 14483321 55305628098268206+0 

. 6411576419708340967533213+0 . 137447407807954792341261 7+0 

.7323364610513124964112125+0 . 13O33O19257384148O0964438+O 

. 8146323700388199325506794+0 .1254359599956809701388411+0 

.8847899667635598184153682+0 .1210294654771558349061820+0 

.94001901 5682225771 1800249+0 .1177316931336763656600726+0 

.9781145826336963259647945+0 . 1155411 781812748548402822+0 

.9975519096742717441995692+0 . 1144492047130534030700801+0 

N=24 

.1058551256599176188873521-2 .1300962489144271767790033+0 

.9300470636577049438343733-2 . 1295533624421600571198363+0 _ 

.2624356686323587846405858'-! .1284714407624273593532948+0 

. 5100890468044952129995145-1 . 12686O974247933OO&423OC0CHO 

.8338425637413211552596118-1 .1247436651916811822325191+0 

. 1228310224356064192750215+0 .1221570801740983595312989+0 I 

.1686928256584271714469156+0 .1191575273509282376358938+0 a m 

. 22CE05 5297233Q0O13O4O5974+O . 1 1 581961 18906842O4759Q616+O ! 

.2765085968459051417058345+0 .1122335336422364381919616+0 ~ a 

.3366578470894900528553694+0 . 1084967032067134824132217+0 s 

.3996397422279969661928756+0 .10470820745803951 50679562+0 ' 

. 4643872630314025802534917+0 . lQ0960789ei316526467Qe278+O 

. 5297973079612462304461 522+0 . 9733577988777187218584927- 1 

. 594749371 5454998136829500*0 . 9390012883373142026738219-1 

.6581251167990389753950110+0 .9070562898232974993311197-1 

. 7188283674669352577185219+0 . 8778983479768393051 TOS4Q1-1 

.7758050151206948194015029+0 . 8517803345876072346931614-1 

.8280623494487953941100798+0 . 8288568870556051525226063-1 

.8746873595551466266667296+0 . 8092032773 72 5685626399568-1 

.9148636O5227877O00D47696D+O .7928680598301284077035461-1 

.9478863111774369315858682+0 . 7798321963200572370647345-1 

.9731753902181926832464217+0 . 7700842802692209837774873-1 

. 9902861 523274990697271967+0 . 7636020069613013542119557-1 

.9989175058794761117551852+0 . 7603662664391345947791117-1 



TABLE 10.- GAUSSIAN QUADRATURE RULE FOR W(T) = /F(l 


TU> 


.2601 0501 65 S99 59 796495 1241+0 
. 7665361 55 163 7932195 540 509K3 


. 1004363127116205741395353*0 
. 3 5630 740941 0391 641 5777329*0 
. 6660620444127490680653213*0 
.9090116995892097449244775*0 


. 5141004963732342471225538-1 
. 194479O457057079315174473+O 
. 3963 5 51 764009910077004924+0 
. 62O631266427OO969402O1025+O 
. 8175026034905643555536350+0 
. 952239148645567 1249552286+0 


.3106074330719916469661529-1 
. 12022010095639767 57871 586+0 
.2559219450409279056324633+0 
. 4207835167132883042571375+0 
. 5942971871195956996833399+0 
. 7556487146290397522942777+0 
. 8860689467717148293265222+0 
. 9706964079107579419636031+0 


. 1483351953900109561287533-1 
. 58429169O6Oe0OD1495 596842-1 
. 1281224330340857762191278+0 
.2196531251903271053302747+0 
.3274580467418496973191181+0 
. 4450681 527883464326806727+0 
. 5655403458716O72212O49566+O 
. 681922789O96176376990O525+O 
. 78756019110361 54371419453+0 
. 87651 53741662579340225392+0 
. 9438262489CE1 344479269577+0 
.9657601749552686610675211+0 


- T 3 ) 


wm 

N=2 

.2433921637378349558374854+0 

.2802066118604639172396218+0 

N=4 

.6122403719459282517444998-1 
. 17758O0214978193857435479+O 
.2Q07456O56425177368489641+O 
. 840491 1 126336892531014518-1 

.2288126739982375588634277-1 
. 7905537440334408253656766-1 
. 1347433688532165788886573+0 
. 1488232296794462160470343+0 
. 1041775419179525725270827+0 
.3391799334451566719142244-1 

Nf=8 

. 1 0628660021 94 5029782891 93- 1 
. 3961 1 53478462413748068179-1 
. 77O1925341757O87955O76093-1 
. 107784452571 582391 5594364+0 
. 1 1 703953987967671 66769444+0 
. 9749236845676136352027743-1 
. 5691187709248358900955109-1 
. 1671108937365476549656322-1 

N=12 

. 3594764872458227443500004-2 
. 138321 1302131 833482397221-1 
.2912099719741149945326390-1 
. 469374704 500360061 1396756-1 
. 63969408056344838724231 18-1 
. 76 58462 591 12467 1 729646734- 1 
. 81384364853788387279781 76-1 
. 76473 77 53 7 59070003 1606439-1 
. 62337O135506784884469287O-1 
. 4219546566973848462657083-1 
.2137283930101251825704453-1 
. 57759373381 58370293288829-2 


T 1 1 ) 


. 865O01476226O4323567O5627-2 
. 342944809*328100048651 1402-1 
. 76O263763258O039O82551238-1 
. 1323678324333478978290352+0 
.2013243878262944039710184+0 
.28O46335160O74173945OO325+O 
. 36701 1 863 1 797044668494902+0 
. 45796629344661 10394845691+0 
. 55O0O429555667O67O6O31493+O 
. 6405934945307417669788932+0 
. 7260939642007900579456573+0 
. 803853435O722832896578696+O 
. 8712944918349106472448407+0 
. 9e6 19260570 59761 272 708501+0 
. 966743391 5380734210899456+0 
.9916172913713338441067073+0 


.3965644162013197670099419-2 
. 1587814766243156590533662-1 
.3548517459811710962022526-1 
. 6248938454626163380634380-1 
.9645340758645334091473429-1 
. 1368270142851174719855317+0 
. 1829565134066113540478940+0 
. 23409624748701 59646926458+0 
. 2894218732448268393856594+0 
. 348O4497492O48409191D6192+O 
. 409028508374321091 739869e+Q 
. 4714026149681493308883371+0 
. 5341804430171462310877473+0 
. 5963737243022884064701209+0 
. 6570D793851O48891Q288533O+O 
. 71 51369459581137177750942+0 
. 76985698306197147867842O1+O 
.8203199096841700892665920+0 
. 8657455757885147535995462+0 
. 905433191 001 8508CE2539853+O 
.9387715194045888297214805+0 
. 9652477945747052088801033+0 
. 9844552040694219635341657+0 
. 9960988368338868164063064+0 


w n 

N=16 

. 1 60424 5208095 5 5 592 7 79H 43 -2 
. 6275572836551 5901 13627559-2 
. 1359823344845988070704896-1 
.2290191525657347519098517-1 
.3328421354236435170087919-1 
.4364046182788257153577264-1 
. 5272484596391097429378456-1 
.5927169031454383441215707-1 
.6218556317536035638097318-1 
.6077226955619545384669695-1 
. 5495057530749846 569776Q2S-1 
. 4537375275123545OO6909O55-1 
.334058972367072769216741 5-1 
.209345035669258545898471 7-1 
. 1004581391673182225368274-1 
. 2629221689261959434535907-2 

U&4 

. 5025641 124260090880707782-3 
. 1969969036556026241 102645-2 
. 440196501008553961 1203784-2 
. 7640091897646723429606940-2 
.1156957254471521079684464-1 
. 16021 53692303 56 501 46032 56-1 
.2079488014255876828804058-1 
. 2565894581451 75CEO49ei094- 1 
. 3035833200238429431232 1 19-1 
.3462148006815756214326463-1 
. 3817423812216867070991985-1 
. 4O7585592370399454533812O-1 
.4215516465965423284615635-1 
. 4220774506303555797326369-1 
. 4084541 121 544222066190239- 1 
. 380998327992535061 S3 52895- 1 
.3411386763245065707276907-1 
.2913940810474475228278918-1 
. 23 523496202667697892 1 9608- 1 
. 1768329733504397974512528-1 
. 1207196573267686197496171-1 
. 7138812400880335464828275-2 
.3287908421919193279712396-2 
. 8397091 1923 58954306721 1 59-3 


NJ 

00 


TABLE 11.- GAUSSIAN 


.2651104697351609990329214+0 
. 93326991 57643 56029 7064340+0 

h*4 

.7757646244400262192063727-1 
.3709225666857466461530439+0 
.7434296492239199226745601+0 
. 9627676763228771239742406+0 

N=6 

.3633643314556056187995606-1 
. 1629224939661174610475239+O 
. 4116135759679648430313137+0 
. 665861954681 1311657248099+0 
. 8787962491335016457236070+0 
. 9922865033909389768639616+0 

*=8 

.2097541112314532654353128-1 
. 1075970682559741279581245+0 
.2515179103879648200009170+0 
. 4328872905305850432124372+0 
. 6254016254744796841728531+0 
. 8001 6032 48197871C332 12378+0 
.9303555157520129386718993+0 
. 9956407313950453451999295+0 

N=12 

. 9562793563395977012124704-2 
.4975672391074898176137756-1 
. 1 1951 1 51 10684994479595657+0 
.2145213878199612456646296+0 
. 32874473450640321 52220645+0 
. 4547364856764484966413293+0 
.5841188562356727415712745+0 
. 7081 58534031 1265674797546+0 
.8183923406800199309194819+0 
. 9O7242486O55345729257OO0O+O 
. 9685676599428793596909655+0 
. 998061 4 56694880761 3 540379+0 


W(I) 


669191302831 7297849469684+0 
1 9528662 5 1 460390CE 551 7871 +1 


1958261 12982136881471 1229+0 
4095325338049126107958979+0 
62273 1 693633691 00031 39233+0 
1393967213871379229883896+1 


. 9E3 7938297440646630582391 - 1 
. 20042 5 7 56 5746 10e 7849663 77+0 
. 290046921064834966851 5900+0 
.3796569304393841575685032+0 
. 5194552986422337267546661+0 
. 1 140093264596650214487619+1 


. 5351492433377449548634005-1 
. 11 893329689581 5293 59O6O0>O 
. 1749650449701 18461 1293713+0 
.2239718463919254186245791+0 
.2733854762616458362964013+0 
. 3369878782605101258880984+0 
.4523102416077387110694063+0 
. 987988845570591468378041 1+0 


. 2447106508548792043843162-1 
. 5563791204311658050665498-1 
. 8415659006599612621617740-1 
. 1095574365039445623726122+0 
. 132321072 4965503908712599+0 
. 1 53636972231 14951 50277439+0 
. 1752678286804981278420637+0 
.1997272841484337614711994+0 
.2311649472792023753151938+0 
. 2784906481818026698393761+0 
. 37050408362 586991 69673861 +0 
.8071217131500678635967205+0 


RULE FOR W(T) = 1/(1 - T 2 ) 3 ' 4 


. 5446744044126310537637219-2 
.2849031076583425111036174-1 
.6910080735458836640883377-1 
. 1258509260932921 175537894+0 
. 1967052422607552814561918+0 
.2790820728668401466076612+0 
. 36993 586391 7976307 1 864689+0 
. 4658648869844356816889626+0 
. 5632397172354233908812175+0 
. 6583455514355559505851456+O 
. 7475307573176861755359536+0 
. 8273 542397052074850064480+0 
. 894724775118061 5424736945+0 
. 9470261 147477487354899174+0 
.9822216061971771322142767+0 
. 996908447 5790109558829880+0 


. 24 50760708120480396760327-2 
. 12869834814090201 10486497-1 
. 3143921672761 707833207207- 1 
. 5786438410895071697799780-1 
. 9171853884313091320475473-1 
. 13245089954272 57308910390+0 
. 1793941044465586669425514+0 
.231 77361 4719e95Q2 59094604+0 
.2887193495618143647956492+0 
. 3492794587895678132322895+0 
. 4124359094813235141O85926+O 
.4771221111344315287763896+0 
. 54224051 1 591 5416945827400+0 
. 6O6682553413579O099O745O8+O 
. 6693477777626274481747269+0 
. 7291633526452684175992687+0 
. 7851030795245740287973647+0 
. 8362055939771172653350511+0 
. 881 59142O73Q90O343534O416+O 
. 920478569674560814381 5727+0 
. 9521963825251747684666929+0 
. 9761973295535414362828386+0 
. 9920600966 537792O880E 5431+0 
. 9995 14421 7484676754187294+0 


U(!) 

N=16 

. 139546276909764108389551 5-1 
.3203063269347150330820596-1 
. 4914687832290443628018975-1 
. 64962 50Qe95O9O999622078 7 5 - 1 
. 794032857 5779 79GE4 704063 5- 1 
.92621 525359X3453732801898-1 
. 1049495283269543467067933+0 
. 1168597544343798691619461+0 
. 1289553223144125961836269+0 
. 14201266944955869796441 1 5+0 
. 1571142009471684659793639+0 
. 17597319O0O61O62845687865+O 
.201 7833074056997940479478+0 
• 241 949571 81 95657666970579+0 
. 321 191 5629010080125732762+0 
.6991489957678901867317743+0 

N=24 

.6284552946755734970665003-2 
. 14533356815702051 10816307-1 
.22 5744891 480300907849e60e-l 
. 3029748073731836034422175-1 
.3763272984763575562688154-1 
. 4454569345249206418404059-1 
. 5105568935384862789699352-1 
. 5713258561 5389521 1 9130666-1 
. 6289204410696678381322444-1 
. 6839085056121 108386553422-1 
. 7372344868471926948617181-1 
. 790005282777 8 277429487019-1 
. 8435009274817467071151954-1 
. 899240051 5018579791477438-1 
.9590525663253643402706934-1 
. 1025294899210576954180014+0 
.1101165264813526279524430+0 
. 1 19132O4O59838848909CB816+O 
. 1303098593856043343533626+0 
. 1449227003493643973814201+0 
. 1654763961418642174583124+0 
. 1979385100553321797914625+0 
.2624375725166264648913186+0 
. 5709713787635603831077756+0 



TABLE 12.- GAUSSIAN QUADRATURE RULE FOR W(T) = 1/(1 - T 2 ) 2 / 3 


Ttl) 


te 2 


W(I) 


T ( I ) 


tel6 


win 


. Z 370t7W4*1903ae7 7204321 4+0 
. 9135412065939290736106419+0 

. te4 

. 7653233930102432946301075-1 
.36561 543277571 1016196021040 
. 7340169991937099640872131 +0 
. 9769204026221 7990787 5 5 51 4*0 

*:6 

.3602196017706051441701710-1 
. 1012 996346320500902 849203+0 
.407600191372651151720031340 
. 6601 50039950238566635330740 
. 0730226249610561 61 96661 OOfO 
. 909562962695339407199576240 


.641984911580965397519369540 
. 146128824640721 599424^1 59+1 


. 1 929856568699660324281308+0 
.3974989635025475272558533+0 
. 568025651 0921998992560396+0 
.9447628865234679328225050+0 


.9155550156709031109846226-1 
. 1978812519807219192116496+0 
. 2823 1 92190986969659268558+0 
. 3564973738548132865082475+0 
.4495745095098818365869159+0 
. 7254453019769770724303975+0 


. 542978841341 69683807328CB-2 
.284006871 5674590759239410-1 
. 6887985272320408385162462- 1 
. 125440955631 7280585382583+0 
. 1960536143168552220171 595+0 
.278147341 75003043 13838CE4+0 
.3686949392533592022938416+0 
.4643196422214430526570766+0 
. 5614230178214642328267213+0 
. 6563228104863825961320366+0 
. 7454005204052867443896808+0 
. 82 52 4466 7 1488822 72901 46 10+0 
. 892787404254 1 401 767806857+0 
. 9454271 736998061 753928660+0 
. 961 131935 1219213389687062+0 
. 9965050085989204636056171+0 


.1391103885278687716824364-1 
.3192635397215304558642838-1 
. 489651 7269825292904576983- 1 

.6465266493602585425443367-1 
. 7885671385605241947004000-1 
.9164837842042515296430830-1 
. 1052553877123435937031 784+0 
.1140165806279100570084840+0 
. 1243557070742238827001678+0 
. 134787681275841 5779436306+0 
. 14597458142460576472929O5+O 
.1588761953193046161651780+0 
. 17513914154939013O1624227+O 
. 1902896087629924391971207+0 
. 238992918935260998554031 1 +0 
.3796250025706120723098011+0 

te24 


te8 


.2084164908665650465654161-1 
. 1068872003788446261851235+0 
.24900020883291 654062 52479+0 
. 42990O5744O1384O7O61863O9+O 
. 62136521950O166971389953O+O 
. 7957227345940077109705212+0 
. 9266520688687292902739341 +0 
. 9940831687467691 740027790+0 


. 5316820947072682296604277-1 
.1179943583965030242324827+0 
. 1726834127852213762564197+0 
. 21 824908385864885820941 52+0 
. 2597310161 548458330358683+0 
.3056748839714324063128131+0 
.3753166298403824934135317+0 
. 60O4555635O842O5773359551+O 


N=12 


. 9522796806409047438131 568-2 
.4954493127894422114136726-1 
.1189901097882763172553211+0 
.2135645322795341638158066+0 
. 327202553223767232753725O+O 
. 4527O486790O99O97673639O8+O 
. 561 0047 1909246O6463397 596+0 
. 7O53320O1245018O094872836+O 
. 81 532071 19649053544225609+0 
. 9046691075507325292131218+0 
. 906724794832572301 3982074+0 
•9973512305128002061930650+0 


. 2436805964941 454829532664-1 
. 55384O93O4310284127512968-1 
. 8366889948233679091726461-1 
. 1085925812746155676184991+0 
. 1303768263966818041 547260+0 
. 1498456496384564155117133+0 
. 1682309488959117823945766+0 
. 1871677878485580625475975+0 
. 20906791 181631 16169261 529+0 
. 2385273341646520210799933+0 
. 2886846337671386774780527+0 
. 45935843201 10012635634962+0 


. 2445710935107822797453948-2 
. 1284318750927283794015242-1 
. 3137358528490702374308480-1 
. 57742290621261 89544 4 048 8 7- 1 
. 91 52265990831601977476978-1 
. 1321645365241536910865975+0 
. 179QD1 87 10891 144679921086+0 
. 231262292631888548167506O+O 
.2880788376684833370054769+0 
.3485037344855510034755054+0 
.4115239971384101771317705+0 
. 47607851 OSCE 120741 683 5926+0 
. 5410762283149457451755300+0 
.6054150878535464422908734+0 
. 6680012531 107167367682701+0 
. 7277682874068713229591969+0 
. 783695888433874 5 1 48GE8214+G 
.8348278303068094393638829+0 
. 8802887808053096284068161+0 
. 9192997031149594162653681+0 
. 951 1915199940648974344861+0 
.9754168352947288898399360+0 
.9915591635768725210049186+0 
. 9993333998143269862485147+0 


.6271585350916622304767417-2 
. 145O088127873926749839482-1 
.2252481411378850588273953-1 
.3022338182712332277208206-1 
.3752258634461899354661572-1 
. 4437812780294266640101332- 1 
. 5077643327871736729338776-1 
. 5673142517374345406468185-1 
.6227988654151781192490377-1 
. 674766137CE9256543O009036-1 
. 7239030182525497726749463-1 
.7710081763657254169110565-1 
. 8169031034946114285342065-1 
. 8628459338064329265368239-1 
. 9097747280793932475322563-1 
.9591955775547927605544475-1 
. 1012929537534091044138457 +0 
. 1O735O10O714333O633277165+O 
. 1144684964495337817743457+0 
. 12327228996654709906541 14+0 
. 1349317518240564311497563+0 
. 152O74O783O4O97O0583O6163+O 
.1828016156571662549162769+0 
. 29Q0e32223486544O96126234+O 



04 

O 


TABLE 13.- ERROR CONSTANTS k N 


Weight 

function 

N 

2 

4 

6 

8 

12 

16 

24 

1 

0.88323-2 

0.35333-4 

0.13925-6 

0.54651-9 

0.83793-14 

0.12817-18 

0.29917-28 

2 

.38568-2 

.15381-4 

.60574-7 

.23766-9 

.36431-14 

.55722-19 

.13006-28 

3 

.83723-2 

.33428-4 

.13195-6 

.51784-9 

.79396-14 

.12145-18 

.28347-28 

4 

.40962-2 

.16231-4 

.63926-7 

.25081-9 

.38488-14 

.58808-19 

.13726-28 

5 

.78509-2 

.31380-4 

.12366-6 

.48528-9 

.74402-14 

.11381-18 

.26564-28 

6 

.43400-2 

.17319-4 

.68214-7 

.26764-9 

.41029-14 

.62755-19 

.14647-28 

7 

.75454-2 

.30175-4 

.11891-6 

.46663-9 

.71543-14 

.10943-18 

.25543-28 

8 

.45125-2 

.18010-4 

.70939-7 

.27834-9 

.42669-14 

.65265-19 

.15233-28 

9 

.18611-1 

.72707-4 

.28401-6 

.11094-8 

.16925-13 

.25831-18 

.60142-28 

10 

.20227-2 

.79014-5 

.30865-7 

.12057-9 

.18397-14 

.28071-19 

.65359-29 

11 

.13386-1 

.51613-4 

.20089-6 

.78342-9 

.11935-13 

.18198-18 

.42341-28 

12 

.11944-1 

.46749-4 

.18286-6 

.71488-9 

.10918-13 

.16668-18 

.38828-28 







TABLE 14.- THREE-TERM RECURRENCE RELATION PARAMETERS FOR WEIGHT FUNCTION 1 
[Zero moment = 0.1402182105325454261175019+1] 


J B ( J) GU> 

1 0.6150128378190335788169563+0 O.* «««••■» o * n m ■ * 00000+0 


2 0.4918^726 7293942056938206+0 0.97199980 72864 766978734869-1 

3 0.4942 420*46 542 7 73 52 7 $6 7 1S2+Q 0. 6^804^960076385^9694835-1 

4 O. 49721 549 73 44C492 68603 5494+0 0. 534716201 7942392756940588-1 

3 0.498335381 73846 7 508 5958096+0 O. 630266558211 7454853634536-1 

6 0.498913830217444 75^3865530+0 0= 62 03C9S2 8 ©91 S3 49954 533 61 5-1 

7 O. 4*55591 580195555707950513+0 O. 62 72 72 76 5633908449263 78 58-1 

8 0. 4S9424 59600267 30460064 59 7+0 0.6266571647154172186701221-1 

9 0.499 5 540 7829 53 5464 3 343 5342+0 0.62^618415588194270540309-1 

10 0.4996442 7681 3 7826402990677+0 0.6259929023366112311948653-1 

11 0.4997096265596906923777585+0 O. <5258016628848901 163648707-1 

12 0.499 7 5846303 1CS5D9334443 183+0 O. 62 566082145291 09093540000-1 

13 0.499793966441 5942918933535+0 0.6255540992892579903320774-1 

14 0.49982 53 5240942270333 12671+0 O. 6254712991341269109622364-1 

15 0.4998488165401011149931981+0 O. ©54057 S8444S89689 178 5706-1 

16 O. 4990678497832O474592O5261+O 0. 6253 5301 €88638101 52S61 730-1 

17 0.4998835016646589778878341+0 0.6253099244959812457538696-1 

18 0.4998965283447 558966728075+0 O. 625274268695787744651 5471-1 

19 0.4999O74858O848277O9476154+O O. 62 S4443 1407201 7 5 53 1043 59-1 

20 0.4999167902564709670214112+0 0.6252192119340475533176225-1 

21 0.49992475*2 5821 298 59329286+0 O. 62 519770410421 ^3630985484-1 

22 0.499931634063911 5 500779 5SOO 0. 6251 792136294655017448687-1 

23 0.499937603 563291 6S9 5228 546+0 0.6251632015092715978563307-1 

24 0.4999428326887550657745987+0 0.6251492432368882826861675-1 

25 0.49994 742 7D48G8 18254084930+0 □ . ^51370024294201042101 103-1 

26 0.4999514889558611242126238+0 0.6251262082814268729665467-1 


TABLE 15.- THREE-TERM RECURRENCE RELATION PARAMETERS FOR WEIGHT FUNCTION 2 
[Zero moment = 0.8413092631952725567050114+0] 


J B ( J) G<J) 

1 0.43959990330 55934134406331+0 *.<■■** ■ *0X00000*0 


2 0.485201 7 i QOl 8394381 5 548000+0 O. 71 1 Ml 5800902764271 103506-1 

3 0. 49 K34Q 72 S 9640061 781 18784+0 0.6442813666042976091603560-1 

4 0.4973594081 7*239831 30 79980+0 0.633296674 5606871531576270-1 

5 Q . 49© 5 1 3 6 962 ©0646300 503901 6+0 O. ^96973 5926022366 53531©34-1 

6 0.499030183939464 782 38 14603+0 0.6260211458494534797182123-1 

7 0, 4 9928 138 5606 56030 58 S 7 402 6+0 0. 6271 C63 54 7685535342098320-1 

8 0.4994 58 5632933 782 3 189S222 5+0 0.6265525449907410750521175-1 

9 0. 499 5 7 73 9^44069926905 71 73+0 0.^6191824^349314627814951-1 

10 0.4996609771701953955917794+0 O. G2 5943 7591 54CSOei 127428059-1 

11 0.4997219907331 51 1660S46697+0 0.6257658533804642161731121-1 

12 0.499767892391 1061 7 S064208+0 O. 62 5633S26 1223 14234 761 7901-1 

13 D.4998CB2923 76 509 57894883 53+0 O. « 5 533388057 5099221001 71 1-1 

14 0.4996311673467357528099563+0 O. ® 54 5 50 12222666 598 59609 52-1 

15 0.4998535090324275498334804+0 O . ^ 5392 7301 5932071 3883S993- 1 

16 0. 499671691 i«O5O0233S92e67+O 0. 62534241 7 5026D2232 13 54261-1 

17 O.4998866058691195975515897+O 0.6253011917047419224366962-1 

1© 0.4998991971771692743253767+0 0. 625266988682871 1531437449-1 

19 0.4999097447421078882977161+0 0. 6252382989 72 52390672266 5 5-1 

20 0.4999187191160587305243560+0 0.6252159979907882585789650-1 

21 0.4999264183447872010118051+0 0. 62 5 1932340 16846 7 50346 163 &-1 

22 0.4999550 730028293368672046+0 0.6251753523559262072777206-1 

23 0. 499938864Q80O47Q&93SGXKK+O 0. 6251 598431194931680803091-1 

24 Q.4S994393463 738931 06544 5 50+0 0.^51463043927422864958548-1 

25 0.4999483994818230691 187482+0 0.6251344159158641393491053-1 

26 0.4999523514007295296044060+0 0.6251239199535311867247599-1 
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TABLE 16.- THREE -TERM RECURRENCE RELATION PARAMETERS FOR WEIGHT FUNCTION 3 
[Zero moment * 0.1311028777146059905232420+1] 

J BtJ) G(J> 


1 

Z 

3 

4 
9 
6 
7 
• 
9 

10 

11 
12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 


0. 9990701 1 7367796103 7199612+0 
0. 498 5 7 5308391 23 1306 594 12 79+0 
0. 4942065962236433906526132+0 
0. 4972508186346309701 777334+0 
0.4983683347062163400825880+0 
0.49892071 1 5762312343521 780+0 
0.4992332802030389296191316+0 
0. 4994272607319624542687284+0 
0.4995558999123650428058766+0 
0.4996455788146774711040877+0 
O. 4997105865356024598802937+0 
0. 49975921 19907668423085779+0 
0. 499796532929108761 5448246+0 
0.4998258013316745562848209+0 
0. 4998491 783012137528924902+0 
0.499868145561 1 728783667530+0 
0.4998837465770967898195050+0 
0.4998967334164776476489372+0 
0. 4999076592322160168824781+0 
0.4999169382228966046379858+0 
0.4999248855154766978370628+0 
O . 499931 7443028827493188126+0 
0. 4999377046879372664677897+0 
0.4999429170091685794121872+0 
0. 4999475014210552731992583+0 
0. 499951 5548854324930980912+0 


o.oo am ' mimiiiiiMHMH O 
0.98061 57552139662646437001-1 
0.6512306687296883760940316-1 
0.6345674644453224600076581-1 
0. 630216622O438324487745683-1 
O. 6282864671 434263341294804- 1 
0. 6272595745632631642853379-1 
0.6266489347407277283216556-1 
0. 6262563614141393223819131-1 
0. 62598907038860555971 701 53-1 
0. 6257988789342640280374204-1 
0.6256587355905226469287555-1 
0. 6255524963080244789721383-1 
0. 6254700407730020482077893- 1 
0. 6254047625922939490837249-1 
0. 6253522OQ26128OD0O4S97432-1 
0. 6253092524613064958977411-1 
0. 6252737090370467398344471-1 
0. 6252439604041 707706777642-1 
□. 62521881 10087462402414351-1 
0. 6251973613204479311178697-1 
0. 6251 7891 7734056964261 5230-1 
0 . 6251 629441264089e741 69945- 1 
0. 6251490182977750387040971-1 
□ . 6251368045618505435808750-1 
0. 62512ea3331d02288933551 57-1 


TABLE 17.- THREE-TERM RECURRENCE RELATION PARAMETERS FOR WEIGHT FUNCTION 4 
[Zero moment = 0.8740191847640399368216132+0] 


j 


6 ( J) 


GCJ) 


1 

Z 

3 

4 

5 

6 

7 

8 
9 

ID 

11 

12 

13 

14 
13 
16 

17 

18 

19 

20 
21 


24 

25 

26 


0. 4493CK 50602 56470777899709+0 
0.4644472950418990698289466+0 
0 . 49 52496639496663 1 3 5967710+0 
0.4975314601266464746424106+0 
O. 4965023969292226020748584+0 
0.4969940496614565716480661+0 
0.4992776422576260706692689+0 
0. 4994561 125571365342526840+0 
0. 4995757074945356397630047+0 
0. 499659760530890784951 7661+0 
0. 49972106666331 70393234334+0 
0.4997672023488136769714523+0 
O. 49960275379476496786781 76+0 
0. 499630738951 3027002553390+0 
0.4998531627246158081290975+0 
0.4996714071878286934799331+0 
0 . 4996664 501 89805331 4 500696+0 
0. 49989890938912099377891 53+0 
0. 4999095771 521 773948082559+0 
0. 4999185758769712546284759+0 
0. 4999262949586761333477443+0 
0.4999329660441485916410501+0 
0. 4999387706320396402101213+0 
0.4999438525599051594134650+0 
0.4999483270082625276689169+0 
0. 4999522870893819907785304+0 


o.oooooQorn~nrrrMTrTTr^ 

0. 7229513301995411333776938-1 
0. 6439933981136391 471 1 74203- 1 
0.6333844246251481963630555-1 
0. 6297404143551978993652384-1 
0.6280415357384917522671346-1 
0.6271182345569966339325942-1 

0.6261968906479339953848994-1 
0. 6259473322475693390116272-1 
0. 6257684674954967965057991-1 
0.6256358957201063222902333-1 
0. 62553490687619471 7321 1291-1 
0.6254562108885519244057762-1 
O. 625393691 5837527579414533-1 
0. 6253432004126267878853528-1 
□. 6253018376844604379688347-1 
0.6252675278900052877356511-1 
0. 6252387536989785377342962-1 
0. 6252143850O9643O1349O9688-1 
0. 6251935661 164539003196196-1 
0. 6251756394597657275221424-1 
0.6251 0009300231 14625219419-1 
0. 6251465232200624994089396-1 
0. 6251346086276909425429012-1 
0. 6251240905454199895408462-1 
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TABLE 18.- THREE-TERM RECURRENCE RELATION PARAMETERS FOR WEIGHT FUNCTION 5 
[Zero moment = 0.1214325323943790805909971+1] 

J BU) C(J) 


1 

2 

3 

4 

5 

6 
7 
a 
9 

10 

11 
12 

13 

14 
13 
16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 


O. 37 73 502691 8962 5764 5091 469+0 
O. 5073244852932461 776984732+0 
O . 4948 79 51 056 76222 733380612+0 
O. 49726720781309814106273000 
O. 4983825039421 185072184619+0 
0 . 4989282968380429702 5121 89+0 
0.4992378482546192918920371+0 
0.4994302188104353046369573+0 
0. 4995579238210856127905567+0 
0. 4996470238734425496061292+0 
0.4997116540121911954372261+0 
O. 4997600227601 1401200e 7351+0 
O. 4997971 63122 782928CE8 1298+0 
O. 4998263008275896318487449+0 
O. 4998495808796588927641945+0 
O. 499868474 7 57 5800e 17293738+0 
0.4998840191953451 561375046+0 
0.499896961 7132075425616461+0 
O . 4999078 52317015 789961 6100+0 
O. 49991 71 0298021 98047639034+0 
O. 49992 502 72262579 7000923 54+0 
0. 4999318670726330257024682+0 
0 . 49993781 1 74 722044830 704 74+0 
0.4999430109281 509055388034+0 
0. 499947584265101 1210207025+0 
0.4999516283295563852513043+0 


O . cn*TwvrrrriTirr^ 

O . 97851 5932049650891 589191 7- 1 
0.66071 51331035211590372786-1 
O . 6342 59 56648239 52 73 5 501 1 40- 1 
0 . 6301 7 73 59 7 7 56086989783538- 1 
0 . 6282606066997644894368888- 1 
O. 627245003942441 4981 486837- 1 
□. 6266390O396645OQ2888O51 71-1 
O. 6262 50274910479483 5022006- 1 
0.6259848117780434682202012-1 
0.6257957836022860350318375-1 
0.62565641 591 59946139222956-1 
O. 6255507132375737305387282-1 
O. 6254686407717625597071 578-1 
0.6254036433358142533188712-1 
0.62 53 51291 4322 70008348 5 744-1 
Q. 6253085044489808305010964-1 
0.6252730860304818867801960-1 
O.625243436O0O6568242235D83-1 
0.6252183662993729371625628-1 
0.6251969796210921368302859-1 
0.6251785882178947886573821-1 
O. 625162657698631 1138886451-1 
0.6251487677616924519194694-1 
O. 6251365841626843977796821-1 
0.6251258384016639071016247-1 


TABLE 19.- THREE-TERM RECURRENCE RELATION PARAMETERS FOR WEIGHT FUNCTION 6 
[Zero moment = 0.9107439929578431044324781+0] 

J B ( J ) G(J) 


1 

2 

3 

4 
3 
6 

7 

8 
9 

10 

11 
12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 


0.4618802153317006116073190+0 
0.4842922474383111053137346+0 
0.4951303458690923596814542+0 
O. 4975060384375699610821943+0 
O. 4984901453737904384305800+0 
O. 4989871213782979101 558020+0 
O. 4992734249441331096310103+0 
0.4994533523708811746413716+0 
O. 4995738038362481431348757+0 
O. 4996583926532631487561469+0 
O. 4997200709624524781621921+0 
0.4997664275917414100886705+0 
O. 4998021494165O9351 5891CE1+0 
0.4998302584402616058791617+0 
0.4996327744167714238694441+0 
O. 4998710689242144546537165+0 
O. 499886186007661 51 58601096+0 
O. 4998987778240754602995538+0 
0. 499909389458421 5274819010+0 
O. 49991841 5487750161 7292742+0 
0.4999261568245507204524516+0 
0. 4999328462309934921 190331+0 
0. 49993866603801 16782913450+0 
O. 4999437607127858799346190+0 
0. 4999482459181 128238046202+0 
O. 49995221 513990677601 1 5839+0 


Q, CT 1 OCPl _ » TTTTTgTTTTTT T' 0 30Q+Q 

O. 7412328435886561499483500-1 
0.6428962495683730122573380-1 
0.6336272305060683443434178-1 
O. 6297744324108638830802408-1 
O.628O65O35143O0O5483O73168-1 
0.627131 58831 12505743978424- 1 
0.6265685484936632039630714-1 
0.6262025937055727058592352-1 
0.6259513511832413969369233-1 
0. 625771405208874561 8903000-1 
0.6256381077933834168882223-1 
0. 62553661 59283641267640380-1 
0.6254575556732714490780611-1 
0.6253947697735688558167979-1 
0.6253440780695731512164820-1 
0.6253025616112264421992180-1 
0. 6252681319936344490558501-1 
0. 6252392630303689392625671-1 
0.6252148184005466529684303-1 
0.6251939379440144269702956-1 
0.6251759608535072589354035-1 
0. 6251 603726861 7026879551 44- 1 
0. 6251467681 103248427978876- 1 
0. 6251348242644722062583706-1 
0.6251242814066083723663851-1 
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TABLE 20.- 


THREE-TERM RECURRENCE RELATION PARAMETERS FOR WEIGHT FUNCTION 7 
[Zero moment = 0.1163592571218269375302518+1] 


j 


B(J> 


G(J) 


t 

2 

3 

4 

5 

6 
7 
6 
9 

lO 

It 

12 

13 

14 

15 

16 
17 
16 

19 

20 
21 
22 

23 

24 

25 

26 


O. 563353679001 4043956734066+0 
0. 511956371 7209448 564 565940*0 
0. 49625404673071 6969604 7499*0 
0.49 72 54306 7 56 51 4 71 19961 466+0 
0.4963921 392124710399060651+0 
0.4969325469450493790249649+0 
0. 4992404466000627510249753+0 
0.4994319027130564147704560+0 
0. 4995590759329695628673962+0 
0.4996476467596316495096792+0 
0.4997122620369225946160946+0 
O.49970O48467564324O3485743+O 
0.4997975222333051213195474+0 
O. 4996265855131419671021 771+0 
0. 4998496103648121979900767+0 
0. 4996686624394869363949409+0 
0.4998841746412091549214587+0 
0. 4996970919025462371016782+0 
0. 4999079624381221807334044+0 
0. 49991 71 969 546 78 730409 7069+0 
0. 499925108062666954075191 7+0 
0.4999319371102442767946124+0 
0. 4999378728269778307367276+0 
0.4999430645148332466406350+0 
0.4999476315359132306546507+0 
0.4999516702393120999680997+0 


0.96845969387166O643O115548-1 
0.6695734049367672681943645-1 
0. 6346714493388579894323679-1 
0. 6301213826076785229810768-1 
0. 6282480423779797464737065-1 
0. 6272365704741117522298901-1 
0. 62663461 79791909180770479-1 
0. 6262468099577302312091488-1 
0.6259823871111705993375620-1 
0. 6257940210634034684810536-1 
0.6256550944897801616402807-1 
0.6255496972716889250477878-1 
□. 6254678429183523106012187-1 
0.6254030053675376051334807-1 
0. 6253507733264366281679306-1 
0. 6253080779626968600244540-1 
0. 6252727307731272196139451-1 
O. 62 5243 1369830966108100339-1 
0.6252181122015177250929386-1 
0.6251967618970646698117800-1 
0.6251784002433503924641996-1 
□ . 6251 62494291 182 53 523981 75- 1 
0.6251486248200589602834210-1 
0. 6251364584069531914566666-1 
0 . 62 512 5727183 53 7021 41 73250- 1 


TABLE 21. - THREE-TERM RECURRENCE RELATION PARAMETERS FOR WEIGHT FUNCTION 8 
[Zero moment = 0.9308740569746155002420145+0] 


J B ( J ) 

1 0.469461 5656345036630611740+0 

2 0.4849526950135912138593998+0 

3 0,4948771 797220851 551234072+0 

4 0.4975175742671236633519254+0 

5 0.4984803228520528929148041+0 

6 0.4989833121879826349802523+0 

7 0. 4992709831252707861053130+0 

8 0. 499451 7627381845929661390+0 

9 0.4995727067142081447920463+0 

10 O.4996576O63O8662990O49961O+O 

11 0. 4997194874448137103223688+0 

12 0.4997659827178376178674O56+O 

13 0.4990018025194919015014299+0 

14 0.4996299827351282878935733+0 

15 0.4998525516811834546477417+0 

16 0.4998709064126078856678562+0 

17 0.4998860346688295244540784+0 

18 0. 49989865081 77181079290609+0 

19 O.4999O9281806225736327411 5+0 

20 0.4999183235784341479388249+0 

21 0.4999260776790245171521407+0 

22 0.4999327775910696071722088+0 

23 0 . 4999386061 2362 53398602670+0 

24 0.4999437061054796138777497+0 

25 0. 4999481994 7634745082 5361 3+0 

26 0.4999521739366385386476169+0 


GU> 

o. am oarrMTTi r g ri^^ 

0. 7547266847071244641 551846- 1 
0. 6423O4342864402O998393976- 1 
0. 633688475361884923O0O6122-1 
0. 6296076834367451093252475-1 
0. 6280770657137456874476197-1 
0. 6271394724628407929925583-1 
0.6265734176943588756216890-1 
0. 6262058778261 776837657233-1 
0. 6259536622512632961448388-1 
0. 6257730934246212654360975-1 
O.625639378293811346O159507-1 
0. 6255375959250571905342461-1 
0.6254583274052183406159435-1 
0.6253953883172798174534847-1 
O.625344581432828973170O786-1 
0. 6253029767076025371695791-1 
0. 625268478312671 6D8CS 79348-1 
0.6252395549657168282015660-1 
0. 62521 50667689005253806993-1 
0. 6251941 510009904994294668-1 
0.6251761449B8218O11492O924-1 
0.6251605329055893731091556-1 
0. 6251469083831203642491290-1 
D . 62 51 3494 77690421 812 58 5691- 1 
0. 6251243907146491366938626-1 
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TABLE 22.- THREE-TERM RECURRENCE RELATION PARAMETERS FOR WEIGHT FUNCTION 9 
[Zero moment = 0.2428650647887581611819942+1] 


J 


B< J) 


G<J> 


1 

2 

3 

4 
3 
6 
7 
a 
9 

10 

11 
12 

13 

14 
13 
16 
17 

15 

19 

20 
21 
22 

23 

24 
23 
26 


431 1849265382984224922325+0 
5344112677183342651850937+0 
4997081663357153245114400+0 
499997992323652642021 6647+0 
5QcmX)753 5 594648049484 54+0 
SOOOOOOOI 321 33 121 42 794695+0 
4 9 99999999 586801654 786577+0 
499999999996862426051 5453+0 
501 llllll U3 JU 2 7 7664 504 764+0 



0011640+0 


4 90999099009 9 9 99999999 9 13+0 

4 90000000090999999999 9968+0 


5GCOnrrici 


0+0 



5r* * ■ i ■ j rrrinrDCOOOOClJUOOO+O 
yii x i'0"*Tnrr)rxTXiioii30000D+0 


121 5894930360789821822357+0 
6302335161643522600332703-1 
6250797750313006342035336-1 
62499961027982358O5729757-1 
624999966221 6613204448769- 1 
62 50000000 7042 7785283 6404-1 
625OO00D00232Q57766949852-1 



624999 9 9 09990 9 99 999997285-1 
hi 



625nrnu rnm 11 1 jl juuooouo- 1 


TABLE 23.- THREE-TERM RECURRENCE RELATION PARAMETERS FOR WEIGHT FUNCTION 10 
[Zero moment = 0.5235987755982988730771072+0] 


j 


B (J) 


G ( J ) 


1 

2 

3 

4 

5 

6 
7 
a 
9 

10 

11 
12 

13 

14 

15 

16 
17 
IS 

19 

20 
21 
22 

23 

24 

25 

26 


O. 5346805956594948912557218+0 
0.4996413756642561247934532+0 
0.4999755312512212612701437+0 
O. 5000000431490767513677482+0 
O. 5000000242857029786729298+0 
O. 5 0 00000000653895318792445+0 
O. 4999999999739740331698369+0 
O.4999999999990O37O1432O393+O 
O. 

O. 

0 . 4 99 9900990099999684488097+0 
0.4 9009000 9 99 909099 9 41 58492+0 
0.5C* 

O.snr 

O. 

O. 

O. 5nrrmTrnriii^ 

O. y m lOTTrimi t lXUi.'iUtW 
O. ^ ■ ■ rinnnfT rir r TnnrTr rT TTT^ 
O. Sn rrrrrrg rrMTCTOOOOOOOOCQ+O 
O. 5ni i rnrnrrrrn u 
0. 5 n r rrrnirmnr n iT nnmrnfn 
O. 5n r im ■ rTrrrrii 1 1 » f. iTJ OLIxao+0 
O. y i iTirrriTr»T« TrmTnTw o 
O. y Tr * ■ i mTiTTm I iimiiBCHO 

o. 



rrxrrrr 

61 7813282621 7020632362457- 1 
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6249999916983040738478722-1 

62 50000000099616865540936- 1 
6240 00 0 0 00007963060999185-1 
624 9 000000009876211407871-1 
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TABLE 24.- THREE-TERM RECURRENCE RELATION PARAMETERS FOR WEIGHT FUNCTION 11 
[Zero moment = 0.2622057554292119810464840+1] 


J 


B'J) 


GCJ) 


1 

2 

3 

4 
3 
6 
7 

e 

9 

10 

11 
12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 


O. 7627597635018131880623260+0 
0.4356406220177258406970294+0 
O. 4836622253812257957504204+0 
O. 4926359637757834902129061+0 
O. 4958160622403241212889954+0 
O . 49 730e 5 7 53 683422 1 98 594969+0 
0.4981164296636529946910406+0 
O. 4986102355901269385187147+0 
O. 4989323521912857650120867+0 
0.4991 541052638607539612490+0 
0.4993132670282797887930142+0 
0.4994313681629815576863545+0 
O. 4995214144758561080843496+0 
0.4995916401095205920752087+0 
O. 49964746473871 566551 16376+0 
0.4996925741142250124147184+0 
0. 4997295457223993006853085+0 
0. 4997602259978321357255499+0 
0.4997859657879120033126338+0 
0. 499807771843331 5107473383+0 
0.4998264069016278370299617+0 
0 . 4998424 573 786985573892202+0 
0.4998563801202517857375893+0 
0. 4998685353030452627082922+0 
0. 4998792099848722599216561+0 
0.4998886352223224496996331+0 


O. OnnrriTTTTTTTTTTTTTTTTTZlXH-O 
O.848642O984932468O6214610O-1 
0.6015784166190404885710378-1 
0.6193678264971025932719100-1 
O . 622 5222 58 503 76234928 122 74- 1 
0. 6236205886612335061 71 5978-1 
O. 624126038343964845673581 7-1 
0. 6243986224336240279938661-1 
0.6245617352484519608591106-1 
0. 624666823492O447603665392-1 
O. 6247383767893202305716068-1 
0.6247892354516655206672295-1 
O. 6248266485212445221838946-1 
0.6248549552240358527784499-1 
0.6248768795508785713376816-1 
0 . 624894200334 1 689 749 5 5 73 42- 1 
0. 6249081 184928942431 655294- 1 
0.6249194681785359645619747-1 
0 . 62492884321 769 76983229232- 1 
O. 6249366755340057480597428-1 
0.6249432852952401654116095-1 
0. 6249489138525996160567481-1 
O. 6249537458977838778488676-1 
0. 62495792468 7D493 52 798 7046-1 
0. 624961 5627029216380232851-1 
O. 6249647492504550799883447-1 


TABLE 25.- THREE-TERM RECURRENCE RELATION PARAMETERS FOR WEIGHT FUNCTION 12 
[Zero moment = 0.2103273157988181391762529+1] 


J B ( J ) GCJ) 

1 a.7131741278126590550076291+0 o . (TixminnnnnnnrrniTiTinw) 

2 0.4574649756003002978073342+0 0.9138266341865190505355904-1 

3 0.4880003168683295789202167+0 0.6214157546585639457435741-1 

4 0.4944457136921354886874163+0 0.6256266089751222097592868-1 

5 0.4968059789822665620923860+0 0.6256299807035465916956906-1 

6 0.4979263010407050897093034+0 0. 6254803891898191389681332-1 

7 0.49854 52766986458567O19637+O 0. 6253641066555803461296587-1 

8 0.4989231367944321599964818+O 0. 625282101 5902222334908550-1 

9 0.4991707046055156317255353+0 0.6252238768292636176628102-1 

ID 0.4993417030237524201212730+0 0.6251815423805393280826464-1 

11 0.4994647555958287820432594+0 0.6251499729997178538489978-1 

12 0.4995562544624473435839419+0 0.6251258760325053868405053-1 

13 0.4996261382321980163249037+0 0. 62 5 107098 75903630600991 49-1 

14 0.4996807182051692635412360+0 0.6250921991884649904453897-1 

15 0.4997241588988035784313898+0 O.625O8O1871899O1279532610O-1 

16 0.4997592984903970030365944+0 0.6250703668503265348254934-1 

17 0.4997881253624322011219338+0 0. 6250622384898620495766689-1 

18 0.4998120661803502762728449+0 0.6250554363219312863216430-1 

19 0.4998321661305429612525275+0 0.6250496878593744372225568-1 

20 0.4998492051120634935611432+0 0.6250447868911064312231361-1 

21 0.4998637746348907397280967+0 0. 6250405751193489093328913-1 

22 0.4998763299112804176901203+0 0.6250369294431442926087396-1 

23 0.4998872258830428922387086+0 0.6250337529972380990634016-1 

24 0.4998967426074243971984465+0 0.6250309687354327723032198-1 

25 0.4999051034483297284078400+0 0.6250285147678126826373071-1 

26 0.4999124883125537396555600+0 0.6250263409262308183342118-1 
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